REPORT  DOCUMENTATION  PAGE 

la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

IK.  RESTRICTIVE  MARKINGS 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

3.  DISTRIBUTION  /  AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
is  unlimited. 

2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 

4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

NSWC  TR  86-328 

5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

6a.  NAME  OF  PERFORMING  ORGANIZATION 

Naval  Surface  Weapons  Center 

6b.  OFFICE  SYMBOL 
(If  applicable) 

Code  R14 

7a.  NAME  OF  MONITORING  ORGANIZATION 

6c.  ADDRESS  (City,  State,  and  ZIP  Code) 

10901  New  Hampshire  Avenue 

Silver  Spring,  MD  20903-5000 

7b.  ADDRESS  (City,  State,  and  ZIP  Code) 

8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

8b.  OFFICE  SYMBOL 
(If  applicable) 

9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

1  8c.  ADDRESS  (City,  State,  and  ZIP  Code)  1 

10.  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM  I  PROJECT 
ELEMENT  NO.  I  NO. 


62633N 


F33327 


11.  TITLE  (Include  Security  Classification) 

Rigid,  Perfectly  Plastic  Analysis  of  Ring-Stiffened  Shells  Under  Dynamic  Loading 


12.  PERSONAL  AUTHOR(S) 

Minos  Moussouros 


13a.  TYPE  OF  REPORT 

13b.  TIME  COVERED 

FROM  TO 

14.  DATE  OF  REPORT  (Year,  Month,  Day) 

1986  August 

15.  PAGE  COUNT 

250 

19  AR3PRACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

This  report  treats  the  case  of  a  long  ring-stiffened  circular  cylindrical  shell  made  of 
rigid-perfectly  plastic  material.  The'JWapproximatl*uJTresca  square  is  used  as  yield  surface 
with  the  corresponding  associated  flow  rule.  The  shell  is  subjected  to  an  external  axisym- 
metrrc  dynamic  overpressure  of  rectangular  form.  Maximum  permanent  deformations  are 
obtained. 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT  21.  ABSTRACT  SECURITY  CLASSIFICATION 

□  UNCLASSIFIED/UNUMITED  Q  SAME  AS  RPT.  □  OTIC  USERS  UNCLASSIFIED 
22a  NAME  OF  RESPONSIBLE  INDIVIDUAL  22b  TELEPHONE  (Include  Area  Code) 

Minos  Moussouros _ (202)394-1681 


DO  FORM  1 473, 84  MAR  83  APR  edition  may  be  used  until  exhausted.  SECURITY  CLASSI 

All  other  editions  are  obsolete.  " 

r^3 _ -  U 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 


NSWC  TR  86-328 


FOREWORD 


This  report  presents  the  first  phase  of  an  ongoing  program  to  establish  the 
permanent  deflections  attained  by  a  circular  cylindrical  ring-stiffened  shell, 
when  subjected  to  dynamic  overpressure. 

The  study  employs  the  ideas  developed  by  Professor  P.  6.  Hodge  30  years 
ago.  It  will  serve  as  a  starting  platform  for  more  complicated  analysis  in  the 
near  future. 

This  work  was  sponsored  by  the  Office  of  Naval  Technology  through  the  Naval 
Surface  Warfare  Center's  Block  Program,  "Explosives  and  Undersea  Warheads” 

(D.  E.  Phillips). 
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INTRODUCTION 


\ 

"ibi 


lis  report  is  the  first  step  towards  an  ongoing  effort  to  obtain  a  closed 
form  or  a  relatively  simple  numerical  solution  to  the  problem  of  a  ring-stiff¬ 
ened,  circular  cylindrical  shell  subjected  to  a  dynamic  pressure  load.  Maximum 
permanent  deflections  after  load  removal  are  obtained. 

The  analytical  solution  of  the  dynamic  problem,  even  under  substantial 
simplifying  assumptions,  can  be  very  complex  and  cumbersome.  In  this  report,  we 
rederive  and  extend  the  analytical  solutions  obtained  by  Hodge — Consider¬ 
able  simplification  can  be  achieved  if  we  are  content  with  numerical  treatment 
of  the  governing  equations  on  the  computer;  fi  ■  ^  . 

This  study  employs  the  assumptions  listed  below. 


1.  (Geometrical  and  Boundary  Conditions  Assumptions)  Rotational  and  axial 
symmetries  are  assumed.  For  ring-stiffened  shells  only  a  typical  half-frame 

spacing (L-j/2)  is  considered.  (This  assumption  implies  a  relatively  long  shell, 
away  from  supports  which  may  affect  behavior.) 


2.  (Loading  Assumptions)  The  dynamic  loading  is  axisymmetric  in  the  form 
of  a  step  function  applied  as  a  pressure  over  a  finite  length  of  time  and 
then  removed  (i.e.,  of  rectangular  type).  The  duration  of  the  load  is  very  short. 
Axial  compression  is  not  accounted  for. 
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3.  (Initial  Conditions  Assumptions)  The  structure  starts  from  rest  with  zero 
initial  velocity  and  displacements. 

4.  (Material  Assumptions)  Material  is  rigid  perfectly  plastic.  The  elastic 
strain  energy  is  much  smaller  than  the  energy  dissipated  through  plastic  defor¬ 
mation.  The  simplest  form  of  the  Tresca  yield  surface  is  employed,  and  no 
residual  strains  are  considered.  No  strain  rate  effects  are  included. 

5.  (Geometrical  Nonlinearity  Assumptions)  There  is  no  geometrical  nonlinearity. 

Small  deformations,  strains,  and  rotations  are  considered.  Transverse  shear 

deformation  is  neglected.  The  employed  strain-displacement  equations  are  linear, 

and  the  equilibrium  equations  are  based  on  the  undeformed  configuration. 

1  2 

The  method  is  not  new;  it  was  first  used  by  P.  Hodge,  Jr.  ’  30  years  ago. 

What  is  new  is  that  we  have  been  able  to  rework  most  of  its  details,  point  out 

any  differences  with  Hodge’s  results,  and  build  the  foundations  on  which  more 

complicated  analyses,  juch  as  by  Jones  (Reference  3)  and  Duszek  (Reference  4), 

1  2 

can  be  completed.  Work  similar  to  Hodge's  original  work  ’  was  also  presented 

by  Kuzin  and  Shapiro, and  Sankaranarayanan^  and  Hodge. ^  Unlike  the  work 
8  9 

by  Hodge  and  Sankaranarayanan,  in  which  the  static  collapse  load  was 
established  under  rigid  plastic  conditions,  the  work  by  Klement^  was  an 
elastoplastic  treatment  of  the  static  collapse  load,  while  References  1  through 
7  were  concerned  with  the  dynamic  case.  Onat,^  Duszek/*’^  and  Lance1^ 
present  analyses  of  load-displacement  predictions  for  post-yield  behavior  under 
static  loads,  accounting  for  changes  in  geometry.  These  analyses  were  based  on 
moderately  large  displacement  theory  and  assumed  rigid  perfectly  plastic 
material.  Furthermore,  Reference  13  presented  a  bounding  principle  for  finite 
deflections  and  static  loadings,  and  compared  it  with  a  solution  by  the  "rate 
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formulation,"  first  proposed  by  Onat  (Reference  11).  Finally,  in  all  fairness, 

14 

the  work  by  Lintholm  and  Bessey  must  be  mentioned,  where  clamped  and  axially 
restrained  beams  were  loaded  impulsively,  and  it  was  concluded  that  the  rigid 
perfectly  plastic  model  was  inadequate  due  to  the  strong  influence  of  elastic 
effects. 

This  report  is  organized  in  the  following  way.  After  a  brief  introduction 
there  is  a  section  of  problem  statement  with  its  proposed  solution  technique. 

The  actual  analysis  pertaining  to  the  obtained  equations  shows  up  in  Appendices  A 
(preliminary  analysis)  through  E.  The  analytical  results  of  this  analysis  are 
summarized  in  Tables  1  through  19  for  the  convenience  of  the  reader.  This  is 
followed  by  a  section  on  results  for  five  cylindrical  shells  and  a  brief  discus¬ 
sion  of  these  computations.  There  is  a  Nomenclature  section  defining  the  terms 
used  and  a  relevant  list  of  references. 
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PROBLEM  STATEMENT  AND  SOLUTION  METHOD 

Consider  a  ring-stiffened,  circular  cylindrical  shell  (Figure  1)  with  con¬ 
stant  frame  spacing  (L^,).  The  shell  has  multiple  frame  spaces.  However,  only 
a  typical  spacing  will  be  analyzed.  The  body  as  well  as  the  externally  applied 
dynamic  loading  are  assumed  fully  axisymmetric  (Figure  2).  The  loading  is  a 
constant  overpressure  of  magnitude  P  and  form  P  ■  P(t)  acting  over  a  time  inter¬ 
val  0  _<  t  <  tg  (Figure  3).  The  time  tQ  is  relatively  short. 

The  material  is  assumed  "rigid-perfectly  plastic,”  [(Figure  A)  References  15 
through  21].  Therefore  there  are  no  deformations  up  to  a  critical  stress.  When 
the  critical  stress,  pressure,  or  load  is  reached,  however,  there  will  be  unre¬ 
stricted  plastic  flow.  In  actuality,  strains  and  deformations  will  not  Increase 
without  bounds  because  buckling  or  fracture  of  the  material  will  instead  take 
place.  In  the  dynamic  case  accelerated  motion  will  be  resisted  by  the  inertia  of 
the  body.  The  rigid-perfectly  plastic  assumption  Is  reasonable,  where  elastic 
deformations  are  very  small  compared  to  plastic  ones. 

To  solve  the  problem,  it  is  required  that  we  employ  some  of  the  ideas  from 
7  1. 5  2  A  25 

plasticity  theory.  ’  ’  ’  Therefore,  we  briefly  mention  that,  when  a  metal 

goes  beyond  its  yield  point,  the  point  representing  the  state  of  stress  in  stress 
space  lies  on  a  surface.  Typicaxly,  we  have  the  Von  Mises  ellipse,  the  Tresca 
hexagon,  and  the  simplification  of  the  Tresca  hexagon,  or  the  Tresca  square. 
Through  a  change  of  variables,  the  stress  space  can  be  written  in  terms  of 
moment  and  force  resultants.  When  more  plastic  flow  occurs,  the  direction 


in  which  this  happens  is  given  by  the  flow  rule.  The  corresponding  strain 
rate  vector  (in  strain  space)  is  such  that  it  is  normal  to  the  yield  surface 
and  is  directed  outwards.  The  employed  yield  surface  (References  7,  15, 
and  19)  is  the  simplified  "approximate"  Tresca  square.  The  associated 
flow  rule  and  normality  conditions  (References  7,  15,  and  19)  are  also 
used.  The  equations  of  equilibrium  (References  1,  2,  7,  and  22)  are  obtained 
in  the  undeformed  configuration.  The  strain-displacement  expressions 
(Reference  22)  used  employ  only  linear  terms.  Therefore,  the  theory  is 
infinitesimal  and  the  analysis  is  simplified  considerably. 

The  boundary,  initial,  and  jump  conditions  (Reference  23)  are  discussed  in 
the  actual  solution. 

The  method  of  solution  is  as  follows: 

1.  Assume  a  "kinematically  admissible”  velocity  profile.  That  is  a  velo¬ 
city  (or  displacement)  distribution  that  satisfies  velocity  (or  displacement) 
constraints;  to  ensure  material  stability  the  total  external  work  by  the  loads 
on  these  displacements  is  positive.* 

2.  Assume  the  portion  of  the  yield  surface  on  which  we  are  on.  We  empx  >y 
”an  associated  flow  rule."  The  "outwards"  normality  condition  of  the  strain  rate 
vector  (written  in  terms  of  displacement  rates)  on  the  yield  surface  leads  to 
certain  conditions  being  satisfied.  These  make  the  mechanism  of  instantaneous 
motion  admissible. 

3.  Satisfy  initial,  boundary,  and  "jump"  conditions  (References  7,  23,  27, 
and  28).  These  last  conditions  are  at  "hinge  circles"  for  shells.  (In  the  case 
of  beams  we  have  just  hinges.) 

A.  Satisfy  the  equations  of  equilibrium  in  the  undeformed  state. 

♦Drucker’s  postulate  on  material  stability.15 
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5.  Verify  that  the  obtained  stress  profiles  (moment  and  force  distributions) 
do  not  violate  the  yield  surface  at  any  point.  If  they  do,  then  another 
"kinematically  admissible"  velocity  profile  must  be  assumed,  and  steps  1  through  4 
repeated . 

The  maximum  "residual”  deflections  are  not  bounds*  of  any  kind.  They  can, 
however,  provide  useful  information  on  damage  (deflections,  not  strains)  by 
dynamic  overpressures  and,  therefore,  serve  as  a  guide  for  the  estimation  of 
deformations . 


*See  upper  and  lower  bound  theorems^  for  static  analysis.  Because  our  problem 
is  dynamic,  estimates  of  this  kind  do  not  produce  bounds. 
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RESULTS  AND  DISCUSSION 


We  summarize  our  analytical  results  in  Tables  1  through  19.  All 

non-dimensional  quantities  that  appear  in  all  the  Tables  (1  through  25)  are 

defined  in  the  Nomenclature  as  well  as  the  Appendices.  Tables  1  through  19  have 

been  organized  with  respect  to  two  non-dimensional  parameters.  The  first  one  is 

a  pressure  parameter,  controlling  whether  the  loading  is  low  or  high.  The  second 

2 

one  is  a  non-dimensional  parameter  c  that  combines  shell  radius  (a),  shell 
thickness  (h)  and  half  frame  spacing  (L).  The  loading  is  termed  "low”  if 

l  +  ^2±P±1  +  ^2  and  "high"  lf  p  >  1  +  ^2  * 

c  c  c 


The  distinction  of  whether  a  shell  is  termed  as  either  short  or  long  depends 

2  ^T  2L^ 

on  whether  the  non-dimensional  constant  c  ■  is  less  or  greater  than  6. 

Tables  1  and  2  pertain  to  short  shells  and  low  loading,  tables  3  through  5  to 
long  shells  and  low  loading.  Tables  6  through  12  to  short  shells  and  high 
loading,  and  Tables  13  through  19  to  long  shells  and  high  loading. 


This  method  was  applied  to  five  cases  of  ring-stiffened  shells,  not  neces¬ 
sarily  representative  of  marine  structures.  Table  20  gives  the  geometrical  and 
material  characteristics  of  these  models.  Model  No.  1  was  subjected  to  an  over¬ 
pressure  of  600  ps 1 •  This  model  was  analyzed  by  using  the  results  of  the  case 
of  short  shells  and  low  loading.  The  applied  overpressure  exceeded  the  critical 
overpressure  (524.569  psi).  Therefore,  plastic  flow  took  place.  Table  21  gives 
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results  for  this  situation.  The  final  permanent  deformations  are  small  In  compar¬ 
ison  to  the  shell  thickness.  At  a  non-dimensional  half  length  of  1.0  (i.e.,  at 

_2 

middle  between  two  stiffeners)  residual  displacement  is  only  0.645  x  10  inches. 

Table  21  also  gives  the  time  for  the  shell  to  come  to  rest. 

Table  22  gives  results  for  Model  No.  2.  In  this  case,  the  length  between 

stiffeners  was  increased  from  3.543  to  10.630  in.  The  analysis  falls  in  the 

2 

long  shells,  high  loading  case,  since  c  (*  28.2546)  exceeds  6.  In  this  case, 
the  applied  overpressure  (600  psi)  also  was  in  excess  of  the  critical  value 
(343.089  psi).  The  residual  displacements  are  an  order  of  magnitude  larger  than 
for  Model  No.  1.  In  this  case,  the  deflection  profile  is  shown  both  as  a  func¬ 
tion  of  non-dimensional  time  (vertically  down)  and  as  a  function  of  distance  from 
the  left  support  (middle  point  is  represented  by  1.000).  Table  23  displays  the 
results  of  Model  No.  3,  identical  to  Model  No.  2,  but  with  lower  applied  over¬ 
pressure,  P  ■  525  psi.  The  residual  displacements  are  smaller  than  for  Model 
No.  2. 

Table  24  gives  results  for  Model  No.  4,  which  differs  in  geometry  from 
Models  1,  2  and  3.  Table  24  shows  that  although  the  applied  overpressure 
(1200  psi)  is  not  much  larger  than  the  critical  overpressure  (986.432  psi),  this 
method  of  analysis  can  lead  to  unrealistic  displacement.  This  is  due  to  the 
omission  of  the  geometrically  nonlinear  terms  in  the  strain-displacement 
relations  (see  Appendix  A,  Equations  (A-12)  through (A-15)]  and  because  the 
equilibrium  equation  (Equation  (A-l)  of  Appendix  A]  was  not  obtained  in  the 
deformed  state. 

Table  25  gives  results  for  Model  No.  5.  In  this  case,  the  overpressure  was 

_2 

reduced  to  987.0  psi  and  its  period  of  duration  increased  from  0.10  x  10  sec. 
to  0.1  x  10  1  sec.  (with  respect  to  Model  No.  4).  We  observe  that  although 
the  exerted  overpressure  (967.0  psi)  is  slightly  larger  than  the  critical  value 
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(986.432  psi),  the  permanent  deformations  obtained  are  larger  than  1  diameter, 
and  there  Is  no  need  to  compute  deformations  to  the  middle  of  the  shell.  The 
same  table  also  gives  the  position  and  the  velocity  of  the  hinge  circle  as  a 
function  of  non-dimensional  time.  Furthermore,  we  observe  that  all  models  come 
to  rest,  as  is  demonstrated  by  the  velocity  of  the  final  point  in  time.  (See 
Tables  21  through  25.) 

In  conclusion,  this  study  has  shown  that  it  is  possible  to  obtain  useful 
expressions  for  permanent  deformations.  However,  we  must  extend  the  method  to, 
at  least,  include  the  effect  of  end  load  [n^  term  in  Equation  (A-30)  which 
will  result  in  a  three-dimensional  Tresca  cube  for  the  approximate  yield 
surface],  and  account  for  geometric  nonllnearltles  in  the  strain-displacement 
relations  to  be  able  to  obtain  useful  expressions  of  residual  deformations. 
Unlike  finite  element  methods,  this  technique  can  provide  very  useful  estimates 
of  residual  deformations  with  minimal  computational  effort,  provided  a  proper 
analysis  has  been  developed. 
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FIGURE  1.  RING-STIFFENED  CIRCULAR  CYLINDRICAL 
SHELL  (FRAME  SPRING  Ly) 


NOTE:  THIS  FIGURE  DISPLAYS  AN  AXISYMMETRIC  STRUCTURE 
AND  LOADING  IN  THE  FORM  OF  INWARD  PRESSURE  P(t) 

FIGURE  2.  CROSS-SECTION  OF  FIGURE  1 
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FIGURE  3.  RECTANGULAR  PRESSURE  DISTRIBUTION 
ACTING  OVER  TIME  tQ 


o 


' - - - 

e 

FIGURE  4.  STRESS-STRAIN  CURVE  FOR  RIGID,  PERFECTLY  PLASTIC  MATERIAL 
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TABLE  1.  SUMMARY,  SHORT  SHELLS,  LOW  LOADING,  1 


SHELL  TYPE 

PRESSURE  LOADING  TYPE 

TYPE 

short  (o  <  C2  <  e) 

CONDITIONS 

0  <  T  <  1 

MOMENT 

RESULTANT 

mx  <T)  “  |j|-  (p  -  1)  -  ijx  3  -  C2  (p  -  1)x2  +  £-£?(p-  1 )  +  3J  x  —  1 

MEMBRANE 

RESULTANT 

=  -1 

DISPLACEMENT 

VELOCITY 

w  { X,  r )  =  |cz  (p  -  1 )  -  2> 

ACCELERATION 

w(x,r)  -  -^g|c2(p-  1>  ~  2JX 

TIME 

o 

N/A 

DISPLACEMENT 
AT  REST 

HAS  NOT  COME  TO  REST  YET 
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TABLE  2.  SUMMARY,  SHORT  SHELLS,  LOW  LOADING,  2 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

■HESBmBH 

CONDITIONS 

1  <  r  <  ~ 

Ko 

MOMENT 

RESULTANT 

mx  (x  ,t  I  =  — 

MEMBRANE 

RESULTANT 

' 

■v“  -1 

DISPLACEMENT 

w  <x,  r)  »  | 

VELOCITY 

w(x,r)=  ~2 (^~  +  ^)T  +  "f"  **j x 

ACCELERATION 

«(x.r)  =  --^(|-+1)x 

TIME  _ 

1  o 

T°  *o  (C2  +  2)P 

DISPLACEMENT 

AT  REST 

13 
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TABLE  3.  SUMMARY,  LONG  SHELLS,  LOW  LOADING,  1 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

long  (c2  >  e) 

LOw(,*i5<P<,.S5) 

CONDITIONS 

0  <  r  <  1 

MOMENT 

RESULTANT 

mx  <x,r  |  =|^{p-1)-lj  x3  —  C2  <p  —  1)  x2  +  [^(p  -  1|  +  3jx  -  1 

MEMBRANE 

RESULTANT 

=  -1 

DISPLACEMENT 

w  ( x, r )  =  _3_|^.(p-1)-ljxT2 

VELOCITY 

w  <x,r>  =  3  |c2(p-1)-2]xt 

2?L  J 

ACCELERATION 

w  { x,  r )  =  3  c2(p-1)-2|x 

2C2  L  J 

TIME,. 

To 

N/A 

DISPLACEMENT 

AT  REST 

HAS  NOT  COME  TO  REST  YET 
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TABLE  4.  SUMMARY,  LONG  SHELLS,  LOW  LOADING,  2 


SHELL  TYPE 


PRESSURE  LOADING  TYPE 


CONDITIONS 


a(r-1)  a 


(C2  -  6) 


1  ^  T  ^  T 

'  *  '  ^  'n 


<u<  i  -vH) 
c 


MOMENT 

RESULTANT 


MEMBRANE 

RESULTANT 


DISPLACEMENT 


VELOCITY 


ACCELERATION 


1  *£(’-£) 


DISPLACEMENT 
AT  REST 


wlx.r  )  °  JL  |£_lp-1)-1 


<  x  <  a  (rc  -  1) 
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TABLE  5.  SUMMARY,  LONG  SHELLS,  LOW  LOADING,  3 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

long  |c2  > 

LOW  (l  +  i_<p  <  1+JL\ 

'  c2  C2' 

CONDITIONS 

(c2  -  6) 

a  ■  v  .  1<t<tu  u  ^  r  <  1 

c 

MOMENT 

RESULTANT 

mK(xrr)  =  1  +  1  r-4x3  +  6(l  +  u)*2  -  12  ux  -  2  (1  -  3uil 

(1  -  u)3L  J 

where  u  *  a  (r  —  1) 

MEMBRANE 

RESULTANT 

1 

ii 

<> 

DISPLACEMENT 

2 

w  (x  ,r  )  = _ -  1)  -  1 

2C2'  2 

-air -1)|  +  11 

a  (1  +a-ar  )J) 

VELOCITY 

ilx.Tl  -  6  /Ml0*-  1  I  l»-») 

q2  \a/  l  6  (1_U)2] 

ACCELERATION 

w (x,  T )  «= 

/  6  \  [1  +a(r-  1)  —  2 x  1  _  , 

'C2'  [l+a-ar]3 

TIME  T 

o 

3( 

ro  -  1  +— 1 

^(p-l)-ll  _ 

2 

c(c  +V6-)  a  '  C  ’ 

DISPLACEMENT 

AT  REST 

W(X,T0)»  3  [-^(P~  1)  ~  l|  X  + 

2C2  l  2  ' 

1  f  6  1 

+«2  U2  1 

IS^HH 
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TABLE  6.  SUMMARY,  SHORT  SHELLS,  HIGH  LOADING,  1 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

CONDITIONS 

'  e*,p-i.  0<r<1  »<«<»„ 

MOMENT 

RESULTANT 

mv  <  x.  t  )  =  2  fiL  -  1 

X  [uo 

3MoB"2[C-6[^]2*6[i]-1 

MEMBRANE 

RESULTANT 

n  .  =  -  1 

DISPLACEMENT 

w(x.r)  =  ~Y  T** 

uo 

: 

VELOCITY 

*  _  <P  -  1)  T  * 
uo 

ACCELERATION 

w-  (P-1)  x 
uo 

TIMEt 

o 

N/A 

DISPLACEMENT 
AT  REST 

HAS  NOT  COME  TO  REST  YET 
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TABLE  7.  SUMMARY,  SHORT  SHELLS,  HIGH  LOADING,  2 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

SHORT  (C2  v<  6j 

CONDITIONS 

2  6  0  <  t  <  1 

U°  =  C2(p-l)  uQ<x<1 

MOMENT 

RESULTANT 

mx  -  1 

MEMBRANE 

RESULTANT 

\  =  -1 

DISPLACEMENT 

w‘y(p"1)t2 

VELOCITY 

W  =  (p  —  1 )  T 

ACCELERATION 

w  =  p-  1 

TIME  r 

o 

N/A 

DISPLACEMENT 
AT  REST 

HAS  NOT  COME  TO  REST  YET 
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TABLE  8.  SUMMARY,  SHORT  SHELLS,  HIGH  LOADING,  3 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

CONDITIONS 

->  6  t  0<x<u  PP 

“  "c2(p->)  ,«=,<,  ' 

MOMENT 

RESULTANT 

*  “j)(t)3  * c2  “2(#  *  f  -  ^){i)  - ' 

MEMBRANE 

RESULTANT 

1 

II 

C 

9  6 

U°  =c2(p -1) 

DISPLACEMENT 


. 


VELOCITY 

w(t)»  (p-r|(*\  ■  C  (p_r  )3/2  x 

W  vT  VT 

ACCELERATION 

•  «rt -£[, ♦£=£•] -4$ 

TfME  r 

o 

AT  TIMEt'  =  ,  TRAVELING  HINGE  u  MOVES  TO  MIDLENGTH  x  =  u  =  1 

p1 

DISPLACEMENT 

AT  REST 

HAS  NOT  COME  TO  REST  YET 
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TABLE  9.  SUMMARY,  SHORT  SHELLS.  HIGH  LOADING.  4 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

SHORT  C2  <6 

C 

HIGH  LOAD  p  >  1  + 

CZ 

CONDITIONS 

u  <  x  <  1 

2  6r  ^  T-  =  —r^ - r  =  - — 

■  -*-■>  Hi-' 

MOMENT 

RESULTANT 

m  =  -  1 

X 

MEMBRANE 

RESULTANT 

r— 

1 

II 

C 

DISPLACEMENT 

w  (r)  =-1-|2pt-  r2-p| 

VELOCITY 

w  (r  )  =  p  -r 

ACCELERATION 

w  (  t)  =  —  1 

1 

TIME  T 

0 

AT  TIME  t'  =  -E_,  TRAVELING  HINGE  u  MOVES  TO  MIDLENGTH  x  =  u  -  1  AND 

Pi 

THIS  REGIME  SHRINKS  TO  ZERO 

DISPLACEMENT 

AT  REST 

HAS  NOT  COME  TO  REST  YET 
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TABLE  10.  SUMMARY.  SHORT  SHELLS.  HIGH  LOADING.  5 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

SHORT  C2  <  6 

HIGH  LOAD  p>  1  + 

C2 

CONDITIONS 

p  >r  t‘  <  r  <  r  r*  -  p 

C2  +  6 

MOMENT 

RESULTANT 

mx  “ 

x3  +  cV  ♦  £_  X  -  1 

MEMBRANE 

RESULTANT 

'V“~1 

DISPLACEMENT 

«-|»*  i 

WHERE  FOR  0  <  x 

«*o<* 

pr.l(c^2)r2j+«i(x,. 

<  uQ  E1  IS  GIVEN  IN  TABLE  11 
<1  E1  IS  GIVEN  IN  TABLE  12 

VELOCITY 

-fi 

ACCELERATION 

TIME  j 

o 

'•-[# 

DISPLACEMENT 
AT  REST 

WHERE  E1  WILL  BE  CALCULATED  EITHER 

T  m  T  BY  TABLE  11  OR  TABLE  12 

0 
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TABLE  12.  En  (x)  FOR  u0  <  x  <  1 


E,  (x)  =  -  — x  p2  Mt  P)c2 

4  (c2  +  6)2 


1  -2p 


C2x  2 

(6 +  02x2) 


+  P 


C4x4 

(eTc  2x2)2 


*#“2« 


+  3C2 

4(6  +  C2) 


3C2  x  ) 
4(6  +  C2x2)  j 


+  x  P 


|  3C2 
((6  +  C2)2 


3C2  x 

(6  +  C2x2)2 
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TABLE  13.  SUMMARY,  LONG  SHELLS,  HIGH  LOADING,  1 


TYPE 

SHELL  TYPE 

AB  PRESSURE  LOADING  TYPE 

LONG  C2  >  6 

HIGH  LOAD  p  >  1  +  — 

C2 

CONDITIONS 

0«t«I 

u  2  -  ® 

c  (P-1)  0  «  x  <  u0 

MOMENT 

RESULTANT 

m„lx,T)  =  2( - 1)  +  1  or  \ 

“°  ) 

=  2(  '  )  -'(  ‘  )  +'(  '  )-’  (  POINTS  ALONG 

°  0  0  V  ar  nw 

MEMBRANE 

RESULTANT 

/  TRESCA 

(  SQUARE 

=  -1  \ 

DISPLACEMENT 

,  ,  1  (P-D* 

w(x,  t)  =  —  - -  t2X 

2  uo 

VELOCITY 

(P-1) 

W  =  -  T  X 

uo 

ACCELERATION 

(P-1) 

W  =  — -  X 

uo 

TIME. 

To 

N/A 

DISPLACEMENT 

AT  REST 

HAS  NOT  COME  TO  REST  YET 
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TABLE  14.  SUMMARY,  LONG  SHELLS,  HIGH  LOADING,  2 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

LONG  C2  >  6 

6 

HIGH  LOAD  p>1+  - 

C2 

CONDITIONS 

..  2  _  6  0  =s  t  1 

^(P-1)  u0<  x  =s  1 

MOMENT 

RESULTANT 

m,  =  1  | 

1 

POINT  B 
i  ON 

MEMBRANE 

RESULTANT 

s 

-6 

II 

>  TRESCA 
f  SQUARE 

1 

DISPLACEMENT 

W  =  1  (P-D  T2 

2 

VELOCITY 

W  =  (p-1)  T 

ACCELERATION 

w  =  p-1 

TIME,. 

T0 

N/A 

DISPLACEMENT 

AT  REST 

HAS  NOT  COME  TO  REST  YET 
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TABLE  15.  SUMMARY,  LONG  SHELLS,  HIGH  LOADING,  3 


TYPE 

SHELL  TYPE 

AB  PRESSURE  LOADING  TYPE 

LONG  C2>  6 

HIGH  p  >  1  +-^T  (p#—  ) 

C2  2 

CONDITIONS 

Oijx^y  I^TigT, 

u  /  6  r  -p  /(p-1)3  ^2  CV  < 

°  V C2(p-1)  1  V  (2p- 3)  j  6  1  j 

©!  defined  in  TABLE  16,  ATTACHMENT  1 

MOMENT 

RESULTANT 

mx  =  -  1 

MEMBRANE 

RESULTANT 

ntp  =  0 

DISPLACEMENT 

1  (P-D  _ 

w  = —  x 

2  u0 

VELOCITY 

o 

ii 

i 

ACCELERATION 

w  =  0 

TIMEX 

To 

N/A 

DISPLACEMENT 

AT  REST 

HAS  NOT  COME  TO  REST  YET 
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TABLE  16.  SUMMARY,  LONG  SHELLS.  HIGH  LOADING,  4 


TYPE 

SHELL  TYPE 

AB  PRESSURE  LOADING  TYPE  1 

LONG  C2>6 

SEE 

HIGH  LOAD  p>  1  +-^—  (p#— ) 

C2  2 

CONDITIONS 

ATTACHMENT  1. 

y<x<u  1<t<t1 

0*-  6,  li- 

C2  [  (p-1)3 

/(p  -  1  )3  \  c2 

T,_p  J 1  (2p-3)  )2  6  6 

(2s/5\-  f 

»*|  0i<min  (l.^5).  if  P>  | 

MOMENT 

RESULTANT 

i  2V3  ^3 

mx(x.r)=  (t;_~j3  [-4x3+6  |u+y)x2  -12  uyx-  C  <  ei  <  ^  ,f  P  <  2 

-  (y + u)(y2-  4yu+  u2)l 

J  •  1 

MEMBRANE 

RESULTANT 

DISPLACEMENT 

SEE  ATTACHMENT  2. 

VELOCITY 

w 

.  ,  (x-y) 

"  lp_T)  0 

ACCELERATION 

M 

•)  1  ,  [u  +  y  -  2x1-  1 

C2/  (u-y3)L  J 

TIME 

To 

N/A 

DISPLACEMENT 
AT  REST 

HAS  NOT  COME  TO  REST  YET 
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TABLE  16.  ATTACHMENT  1 


n/3 


V  =  uo-0+  „ 
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TABLE  16.  ATTACHMENT  2 


V  <x<u. 


1  (P-D  J 

w(x,t)=  —  — —  x+ 


C2  (p-1)3 
6  (2p-  3) 


(  2  \/3  \ 

(x-uo)(0-uo)  +  ~  (02~uo2)  -  (0-uo)  lo8e 

\u°  c  / 

2N/3) 

- 

V  0  c  / 

£  (•♦¥N 


^q+2n/3\|  /..  x  M\ 


:)K- 


+  ~C~) 


i°g0 


-(•-¥)«.  I  (-¥)  |- •(■•-¥)».  1 1--5#)!! 


U.  <  X  <  u 


w(x,r)  =  —  (p-  1)  + 


_ 

1/ 

C2  (p-1)3 

1  ,  V^3 

(X-Uo)(0-Uo)  +  —  (02-Uo2)  -  “20"  (0-uo>  loflB 

C 

O 

1 

O 

6  (2p-  3) 

(„  +  2^\ 

- 

V  0  c  / 

2C  I  \  C 


( +  ^ 

\u° 


c  / 


2^ 


o  c 
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TABLE  17.  SUMMARY,  LONG  SHELLS,  HIGH  LOADING,  5 


SHELL  TYPE 

AB  PRESSURE  LOADING  TYPE 

TYPE 

LONG  Cz>  6 

high  p>i+|t 

CONDITIONS 

u<  x<  1 

MOMENT 

RESULTANT 

mx=i 

MEMBRANE 

RESULTANT 

«V  =  -1 

DISPLACEMENT 

SEE  ATTACHMENT  1. 

VELOCITY 

w=p-r 

ACCELERATION 

w  =  - 1 

TIME  _ 

To 

N/A 

DISPLACEMENT 

AT  REST 

HAS  NOT  COME  TO  REST  YET 

NSWC  TR  86-328 


TABLE  17.  ATTACHMENT  1 


u<x<l  u=u*  at t=t*,  0  =  Q# 
w(x,t)  =  <P“  1)  +  P  [t~t*]  -  [t2-t*21  + 
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TABLE  18.  SUMMARY,  LONG  SHELLS,  HIGH  LOADING,  6 


TYPE 


CONDITIONS 


SHELL  TYPE 


LONG  C2  >  6 


PRESSURE  LOADING  TYPE 


6  3 

HIGH  p>1  +  —  pi  — 


+  V3  . 

yi  =  uo~^  +l£-,09e 


(&1  +  c  )  (u0  -  c  ) 


(0,-  (u  +  ili) 


2V5"  w 

r,  ^  r<r  9,  =  - 


ih  (A)  + 


r  2v3  -i  ■ 

1  +  - coth  (A) 

L  Cu  J 


0  <  x<  y  y  < 


v  o 

j(p-D3  r  c2  l(  % 

r  .=P-< -  2 - e2\  , 

1  |  (2p  —  3)  .  6  J) 


,-sfi 


MOMENT 

RESULTANT 

mx  =  -1 

MEMBRANE 

RESULTANT 

o 

(1 

c 

DISPLACEMENT 

.  .  1  <p-1) 

w(x,t)  = - X 

2  % 

VELOCITY 

w  =  0 

ACCELERATION 

w  =  0 

TIME 

To 

...  .  MM 

*  H.*f] 

DISPLACEMENT 

AT 

REST 

W<X,  ro) 
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TABLE  19.  SUMMARY,  LONG  SHELLS,  HIGH  LOADING.  7 


SHELL  TYPE 


LONG  C2  >  6 


PRESSURE  LOADING  TYPE 


HIGH  p>1  +  — 


2v3 

V3  ,  (6i  +  C  )  (uo_  C  * 

y  =  u„  -0  + - loge - 

2C  ,*  2v^  /  +  2vT, 

(Q - )  (u  +  - ) 

1  r>  O  /■» 


|H~ 

Cv/p^l 


T^  <  r^'( 


2^  LlAlwl  C(1"' 

_  til  L_ _ CugJ  A  -  -—r- 


1  +  - coth  (A) 

Cu_ 


y  <  x<1  y  <  1 


,  6  1 

V=  1 - — 

C2  U  _ 


l(p-1)3  -  c2 

=  p_  { -  2 - 0 

1  |  (2p  —  3)  L  6 

'  1  H')  L 


(i-y,)2J  (p-t,) 


6  1  1/1-' 


(T"ri)  +  Vi 


y(T)  =  h_ 


C<1-Vi  P-r 


=  1+ - ] -  £— 4x3  +  6  (1  +  y)  x2  —  12  yx  —  2  (1  —  3y)] 

(1  -y>3  _  _  , 


nlp  =  -1 


1  (p-1) 

wly,  t)  =  —  -  y 

2  u„ 


(r^D,  C2  A^t-t^-Ib^t2-^) 


1  (Cl"Bir)  (1  -  x)  (  1 

6  -loge  - - -  +-—  7— 


Bi  (^-8,7,)  B, 


~V  C|~Vl 


(SEE  NEXT  PAGE 
FOR  MORE 
INFORMATION.) 


(p-Ti)(1-y1) 


,  ,  u  ,  n/6"  u  ,  .  ^  (SEE  NEXT  PAGE  FOR 

w(yQ  .  to)  when  y  =  y<j  =  i  -  _  the  veloaty  becomes  zero.  MORE  INFORMATION.) 
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TABLE  19.  (Cont.) 


r.  6  1  H-y,) 

ip^ 

±  — !— 1 
L  c2  n-ytfj 


C,=  1-y,+ 


1 

H-y,)2_ 


(i -y,) 

(P-T,) 


(1-y,)  (p—  r ■, ) 

1_  [ca(1-y/-6] 
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TABLE  20.  (CONTINUED) 


TABLE  21.  RESULTS  FOR  MODEL  NO.  1  (SHORT  SHELLS.  LOW  LOADING  CASE) 
(P  =  600  psi) 
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TABLE  23.  RESULTS  FOR  MODEL  NO.  3  (LONG  SHELLS,  HIGH  LOADING  CASE) 
(P  =  525  psi) 
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NONDIMENSIQNAL  HALFLENOTH  RESIDUAL  DISPLACEMENT  NOND I MENS I ONAL  DISPLACEMENT  DISPLACEMENT/THICKNESS 
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NOMENCLATURE 


A 


a 

c,  C 


% 

F(t) 

F(t) 

H 

h 

uj) 


2  o  H  t2 

_ y  ° 

2 

a  s 

Mean  Shell  radius  (In.) 

/2L 

Non-dimensional  parameter  defined  as  - 

✓ah 

Rate  of  plastic  work  done  (plastic  dissipation)  /  unit  area 

(Lb-IN/IN3*SEC) 

Strain  rate  vector  with  components  (e^,  icx) 

Total  axial  strain 

Total  circumferential  (hoop)  strain 

Function  with  respect  to  non-dimensional  time  variable  t 

Derivative  of  F(t)  with  respect  to  t 

Half  thickness  of  shell  (in.) 

Total  thickness  of  shell  (in.) 

Unit  Vectors  on  (n$,  n^)  or  (e^,,  tcx)  space 
(See  Appendix.) 
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U)  n(2) 

SA  ’  ~A 


„(1)  n(2) 

~B  ’  ~B 


NOMENCLATURE  (Cont.) 


■  Bending  curvature  in  axial  plane  (1/in.)  (ic  *  Nondimensional 

•  1  H  .  * 

bending  strain  rate  and  *  —  k^) 


Bending  curvature  is  a  plane  perpendicular  to  the  shell  axis 
(1/in.) 


Yield  stress  in  pure  shear  defined  in  terms  of  Oy  yield  stress 
as  igoy  (psi) 


Half  length  of  shell  (from  stiffener  to  midbay)  (in.) 


Total  length  of  shell  (L<j  ■  2L)  (in.) 
Mass  of  section  analyzed  (Lb-sec^/in.) 


Bending  Moment/length  of  section  of  shell  perpendicular  to 
x-axis  (longitudinal  axis).  (Lb-in. /in.)  It  is  defined  as 

rH 

J  z  oxdz,  where  z  is  dummy  variable  across  thickness 
-H 


Non-dimensional  bending  moment  based  on  M, 


OyH^  (Lb-in. /in.) 


Unit  normals  at  point  A  of  yield  surface  (See  Appendix) 


Unit  normals  at  point  B  of  yield  surface  (See  Appendix) 


axdz  (lb. /in.) 
-H 


2oyH  (lb. /in.) 
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NOMENCLATURE  (Cont.) 

Membrane  force/unit  length  of  axial  section  and  a  section 
perpendicular  to  the  axis  of  cylindrical  shell  (defined  as 

rH 

J  a^dz)  (lb. /in.) 

-H 


Non-dimensional  membrane  force  variable  based  on  N$ 


External  pressure,  function  of  time,  i.e.,  P  -  P  (t)  (pel). 
In  this  report  it  is  taken  as  a  constant 


Non-dimensional  pressure  parameter  I  *  P 


Non-dimensional  collapse  load  parameter  defined  by  p  ■  1  +  — 

°  C2 

Non-dimensional  collapse  load  parameter  defined  by  p^  *  1  +  — 

c 

Surface  density  of  material  (mass/unit  area)  (lb.-sec.2/in.^) 


Time  variable  (sec.) 


Time  period  over  which  pressure  loading  P(t)  is  acting  (sec.) 


Axial  distance  (in.) 


Non-dimensional  variable 


(-?) 


Longitudinal  displacement  (in.) 


Non-dimensional  axial  displacement  ^  jr  j 


Velocity  of  propagation  of  traveling  hinge  (in. /sec.) 


Radial  displacement  (in.) 
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NOMENCLATURE  (Cont.) 


Non-dimensional  radial  displacement  variable  defined  by 


sa 


4<Ht 


W 


sa 


2c  Ht 

y  o 


-  ”  —V 
2  aA 


Derivative  with  respect  to  non-dimensional  axial  variable  x 


Time  derivative  with  respect  to  non-dimensional  time  variable  t 


Square  brackets  indicate  jump  conditions  of  a  variable  (i.e., 
difference  of  its  value  at  two  points,  one  to  the  left  and  the 
other  to  the  right  of  the  point  in  question) 


Arbitrary  positive  constants  used  in  Appendices 


Axial  strain  in  middle  surface  of  skin  (otherwise  known  as  axial 
inplane  strain) 


Circumferential  strain  in  middle  surface  of  skin  (otherwise 
known  as  hoop  strain) 

Non-dimensional  stress  resultant  vector  (n^,  mx) 

Yield  stress  (psi) 


Non-dimensional  time  variable  (t  “  t/to) 
Non-dimensional  time  at  which  shell  comes  to  rest 


Non-dimensional  time 
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This  Appendix  presents  the  derivation  of  the  developed  solutions  in  a 
detailed  way.  The  circular  cylindrical  shell  is  ring-stiffened,  fairly  long, 
and  subjected  to  an  external  dynamic  pressure  load  P,  which,  in  itself,  is 
assumed  axisymmetrlc  and  whose  form  over  time  is  rectangular.  Because  the 
material  is  assumed  as  rigid  perfectly  plastic,  there  will  be  no  deformation 
unless  the  external  load  P  exceeds  the  critical  collapse  load  Pq  (References 
A— 1  and  A— 2 ) . 

For  the  purpose  of  solution  we  must  have  at  our  disposal: 

1.  The  equations  of  equilibrium 

2.  Strain-displacement  relations 

3.  A  yield  condition  and  associated  flow  rule. 

A.  Conditions  of  continuity  of  certain  variables,  and  initial  and  boundary 
conditions. 
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The  analysis  follows  the  five  steps  of  the  problem  statement  listed  in  the 
main  text  for  all  four  cases  listed  below.  (It  will  become  clear  after  each 
solution  in  later  appendices  why  such  division  takes  place.) 

1.  Short  shell  and  low  intensity  loading. 

2.  Long  shell  and  low  intensity  loading. 

3.  Short  shell  and  high  intensity  loading. 

4.  Long  shell  and  high  Intensity  loading. 

EQUATIONS  OF  EQUILIBRIUM  AND  STRAIN  -  DISPLACEMENT  EQUATIONS 


If  we  consider  a  section  of  a  shell  element  of  axial  length  dX  and  cir¬ 
cumferential  length  ad*  (where  a  stands  for  the  shell  radius)  of  mass  M 

(Figures  A-l  and  A-2)  acted  on  by  external  pressure  P,  we  can  obtain  Timoshenko's 

A-3 

equations  of  equilibrium  in  the  undeformed  configuration.  The  external  loads 
for  this  axisymmetric  case  are  resisted  by  internal  bending  moment/unit  length, 

M  ,  and  membrane  force/unit  length,  N  .  The  equilibrium  equation  is: 


M^-rW  -  a  dodx  rM  +  a  d*dX—  +  a  dX  d*P  (A-l) 

3tZ  3X^  X  a 
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Replacing  the  mass  per  unit  area  s  in  Equation  (A-l)  by 


M 

a  d$dX 


(A-2) 


we  have 


s 


a  1 

+  P 

3X2  1  a  * 


CA-3) 


In  the  paragraphs  that  follow,  equations  will  be  written  in  non-dimensional 
form.  The  objective  is  to  obtain  certain  relations  between  certain  groups  of 
parameters,  on  one  hand,  while  on  the  other  (it  is  easier),  to  work  with 
non-dimensional  groups. 

We  introduce  the  following  notation: 


T  -  t/tQ 


(A-4) 


x 


(A-5) 


m 

x 


(A-6) 


A-4 


k 
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\ .  A 

Nc  2ayH 


(A-7) 


P  " 


P_  a_ 

ct  2H 

y 


(A-8) 


w  -  -  .W  (A-9) 

2 a  Ht 

y  o 


2 

c 


( A-10 ) 


where  L  is  the  half  length  of  the  shell,  H  is  the  half  thickness,  ay  is  the 
yield  stress,  t  is  the  duration  (secs)  of  the  rectangular  pressure  loading, 
and  X  is  the  axial  distance. 


Using  these  equations.  Equation  (A-3)  is  transformed  to 

2  2 

— 9  3~om-  +  n*  +  P  "  3— 9W  ”  0  (A-ll) 

2cZ  3x2  *  *  3t2 

This  represents  the  equilibrium  equation  of  an  element  in  the  undeformed  confi¬ 
guration.  It  is  written  in  terms  of  two  generalized  stress  resultants  only, 

m  and  n  .  This  is  due  to  the  axisymmetric  nature  of  the  problem  and  the 
x 

absence  of  axial  loads.  These  resultants  appear  in  the  plastic  dissipation  rate 
D,  depending  on  the  stress-state  and  the  yield  surface  (References  A-l  and  A-4). 
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DISSIPATION  EQUATIONS 


The  purpose  of  the  analysis  that  follows  is  to  obtain  the  plastic 
dissipation  rate  in  terms  of  the  strain  rates  associated  with  the  generalized 
stress  resultants.  The  step  is  vital  because  it  identifies  the  strain  space 
corresponding  to  the  space  of  "stress  resultants"  of  the  yield  surface 
(References  A-l  and  A-4) « 

To  further  the  analysis,  we  state  the  strain-displacement  equations  from 

Reference  A-3.  The  inplane  axial  (e  )  and  circumferential  strains  (e.) 

x 

as  well  as  curvatures  <  and  kx  are: 

x  * 


e 

x 


3U 

3X 


(A-12) 


e 


♦ 


(A-13) 


K 

X 


(A-14) 


K 


♦ 


-  0 


(A-15) 


where  U  and  W  are  the  axial  and  radial  displacements  at  the  midsurface  of  the 
shell,  independent  of  the  thickness,  variable  z.  They  depend  on  X  only  due  to 
axisymmetry  (Reference  A-3). 
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The  total  strains  ex, 


are  given  by 


e 


x 


e 

x 


+  Z  K 


X 


12  +  zl— w 

3X  3X2 


(A-16) 


+  Z  K 


(A-17) 


He  Introduce  a  non-dimensional  thickness  variable  as 


C 


(A-18) 


2 a  Ht  2 

and  A  *  — (A-19) 

as 


Naturally,  A  is  a  positive  quantity. 


We  can  write  Equations  (A-16)  and  (A-17)  as 


3u  ,  „ 

—  +  cHk 
3X  X 


(A-20) 


e.  *  -  Aw 
v 


(A-21) 
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and  Equations  (A-12),  (A-13),  (A-14) ,  and  (A-15)  as 


e 

x 


3u 

3x 


(A-22) 


e 


<t> 


■  -  Aw 


(A-23) 


K 

X 


1_  A 

2  H 
c 


(A-24) 


<  -  0  CA-25) 

To  identify  which  parts  of  the  strain  rates  play  a  part  in  the  analysis,  we 
must  compute  the  energy  dissipated  pJ nstically  (plastic  dissipation  rate  D). 


From  first -principles  (References  A-l  and  A-4),  the  internal  rate  of 
plastic  work  D/unit  area  in  terms  of  stresses  and  strain  rates  Integrated  over 
shell  thickness  2H  is 


A 


-H 


e  +  eA)dz 
xx  $  $ 


(A-26) 


If  we  introduce  the  non-dimensional  parameters  n  ,  m  ,  and  n  by 

XX  9 


means  of 


N  -  N  n  -2c  Hn 


(A-27 ) 
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N. 


N  n. 
o  4> 


2a  H  n 
J  t 


(A-28) 

i 


M  -Mm  -  a  H2m  -  ~  N  m 
x  ox  y  x2ox 


(A-29) 


Equation  (A-26)  can  be  expressed  as 


D  -  e  N  +  k1  M  +  eA  N 

XX  XX  $  <t> 


e  n  +  fc1  m  +  ea 
XX  XX  <1 


N 


(A-30) 


where 


.1 

K 

X 


A 


3 


2 

w 


and  this  is  the  only  time  (numerical)  1  does  not  denote  differentiation. 


In  the  present  case  there  is  no  axial  load  (nx  -  0)  and  Equation  (A-30)  becomes 


.  -1 

e4  *x  ®X 


Nr 


(A-31) 


■  U 


in  terms  of  two  generalized  strain  rates  (e.,  k  )  and  stress  resultants  (n.  ,  m  ). 

<t>  X  U>  X 


We  observe  that  we  can  write  the  generalized  strain  rate  in  terms  of  a 
two-component  vector  je,  in  strain  space,  of  the  form 


i-  <v  s’ 


(A-32) 
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or  in  view  of  Equations  (A-23)  and  (A-24) 


(A-33) 


where  A  is  a  constant  given  by  Equation  (A-19).  Furthermore,  the  stress 
resultants,  which  contribute  to  plastic  dissipation,  form  a  two-component 
stress-resultant  vector,  0  of  the  form 


o  *  (n  ,  m  )  (A-34) 

V  X 

In  the  analysis,  a  yield  surface  and  associated  flow  rule  must  be  used.  The 
yield  surface,  which  satisfies  the  Tresca  yield  criterion,  has  been  developed  by 
Drucker,  Hodge,  and  Onat.  (See  References  A-4  and  A-5  for  more  details.)  For 
the  present  case,  the  exact  condition  will  be  replaced  by  the  simplified  Tresca 
rectangle  ABCD  (Figure  A-3).  This  locus  envelopes  the  exact  curve  from  the 
outside.  According  to  plasticity  theory,  the  strain  rate  vector  e  must  always 
be  normal  to  the  yield  surface  and  directed  outwards  during  plastic  flow,  except 
at  the  corners,  where  e  must  lie  within  the  space  described  by  the  (two  or 
three)  normals  to  the  yield  surface  there.  Furthermore,  the  yield  surface  must 
be  convex.  (See  References  A-4  through  A-10.) 

The  flow  rule  (outward  normality  to  yield  surface),  therefore,  determines 
the  direction  of  growth  of  plastic  flow,  which  takes  place  along  the  gradient  of 
the  yield  surface.  For  perfectly  plastic  materials,  however,  the  magnitude  of 


the  strain  rate  cannot  be  determined. 
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Furthermore,  notice  that  on  a  straight  side,  while  the  strain  rate  vector 

does  not  uniquely  determine  the  state  of  stress  [l.e.,  components  (n  and  m  )] 

$  x 

the  plastic  dissipation  Is  unique.  [See  Hodge  (Reference  A-5),  Chapter  8,  pp. 
195-201.]  Accordingly,  the  strain  rate  vector  e  must  be  of  the  form  Indicated 
In  Table  A-l. 

Here  we  explain  how  Table  A-l  is  constructed. 

1.  Region  AD 

We  consider  the  region  AD  (except  the  two  end  points  A  and  D  which  must  be 
treated  separately  as  they  constitute  "corners"  on  the  yield  surface).  The 
equation  of  the  outward  normal  to  AD  is: 

%  "  C1  -1  i  <=1  £  1 

(A-35) 

Its  slope  is  on  the  (n  ,  m  )  plane.  The  strain  components  of  e  are 

v  x  ~ 

given  by  Equation  (A-33).  Since  mx  *  -1,  this  means  that 


2c 


»w  <0 
2  xx  — 


and  since  A  >  0  ,  this  leads  to  w  >  0. 

xx  — 


(A-36) 


Noting  that  n  is  arbitrary  [to  the  extent  that  it  lies  between  (-1,  1)]  and 
$ 

the  slope  of  the  e  vector  is 


A-ll 
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(A-37) 


this  can  only  be  satisfied  by  w  ■  0. 


2 .  Point  A 


At  point  A  there  are  two  outward  normals. 


a.  Normal  to  AD  ■  nf^  =  (0,  -1) 

“''A 


or  in  (i,  j)  notation 


(A-38) 


b.  Normal  to  AB 


(-1,  0)  -  -i 


As  already  mentioned  (and  discussed  in  References  A-l,  A-4,  and  A-5),  at  a 
corner  of  the  yield  surface  the  strain  rate  vector  e  must  be  a  positive  linear 
combination  (a,  8  >  0)  of  the  unit  exterior  normals,  i.e., 


(-8,  -a) 


(A-39) 
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By  Equation  (A-33),  however,  we  note  that 


k  ’  A(f*’ "  ^2*=)  ‘  ("B’  'a) 


(  A-40 ) 


i.e.,  -B,  -a  <  0  and,  hence, 


-w  <  0 


(A-41) 


2c 


rw  <  0 
2  xx  - 


or 


w  >0 


w  >0 

XX  — 


3.  Point  B 


Similarly,  at  point  B  the  outward  normals  are  given  by 


a.  Normal  to  AB  at  B:  ■  (-1,  0)  ■  -  i 


(A-42) 


b.  Normal  to  BC  at  B:  tu  *  •  (0,  1)  -  j 


42) 


Since  e  must  be  a  positive  linear  combination  of  n^  and  n^,  we  have 


a  n^1)  +  B  n£2)  -  (-a,  8) 


(A-43) 
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Comparing  Equation  (A-43)  with  (A-33)  we  deduce  that 


-w  <  0  or  w  >  0 


(A-44) 


— sw  >0  or  w  <  0 

_  2  xx  —  xx  — 

Zc 


(A-45) 


Table  A-l  can  be  completed  by  similar  arguments. 


CONTINUITY  CONDITIONS 


The  shell  in  question  is  made  of  material  without  voids  or  other  sources 

3  w 

of  discontinuity.  Consequently,  the  radial  displacement  w  and  velocity  are 
continuous.  If  we  denote  "jump”  conditions  through  brackets,  then 


(A-46) 

(A-47) 


In  addition,  the  bending  moment,  m  ,  membrane  force,  n  ,  and  shear  force 

X  y 

distributions  must  be  continuous  across  the  length  of  the  shell. 


Finally,  we  need  to  quote  some  results  from  the  travelling  hinge  theory  of 
A-ll  A-l 2 

Lee  and  Symonds.  *  Recall  that  a  hinge  forms  at  some  point  on  the 

structure,  when  the  limit  moment  has  been  reached.  In  structural  dynamics,  this 
hinge,  hinge  line,  or  hinge  circle  may  either  be  stationary  or  travelling  with  a 
velocity  or  propagation  V.  In  the  present  case  of  our  cylindrical  shell,  note 
that 


A-14 
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a.  If  the  hinge  is  stationary  (V 
discontinuous ,  i .e . , 


0),  then  the  slope  may  be 


'M]*°  «» 


(A-48) 


b.  If  the  hinge  is  travelling  (V  ^  0),  then  the  slope  is  continuous 


(A-49) 


INITIAL  AND  LOADING  CONDITIONS 

The  shell  is  initially  at  rest.  Up  to  time  t*0(x*»0)  there  is  no 

applied  load.  Suddenly,  at  t  -  0  (t  -  0)  a  rectangular  pressure  pulse  of 

magnitude  P  (p)  is  applied  on  the  cylindrical  shell.  The  pressure  pulse  acts 

over  a  time  t  (x  *  1)  and  then  becomes  zero.  Consequently,  the  initial 
o 

velocity  and  displacement  conditions  are 

w(x,  0)  -  0  (A-50) 

w(x,  0)  -  0  (A-51) 

At  some  time  t^  >  t^  (x^  >  1)  the  body  is  brought  to  rest,  i.e., 

w(x,  x^)  *  0  (A-52) 
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The  loading  conditions  are 


p(x,  0)  -  0 

p(x,  t)  -  p 

for 

0  <  T  <  i 

(A-53 ) 

p( x,  t)  -  0 

for 

T  >  1 

where  p  is  a  constant. 


BOUNDARY  CONDITIONS 

With  reference  to  Figure  A-l,  the  first  ring  stiffener  is  located  at  x  “  0 
(X  *  0),  and  the  second  one  (end  of  bay)  at  x  *  2  (X  ■  2L).  At  the  first 
stiffener  (x  *  0)  the  radial  displacement  and  velocities  are  zero  for  all  times 
t  ^  0 1  i.e. , 


w(0,  t)  -  0  (A-54) 

w(0,  t)  -  0  (A-55) 


From  the  discussion  in  "conditions  of  continuity,"  we  conclude  that  whether  the 
slope  of  a  hinge  is  continuous  or  not  will  depend  on  whether  the  initially 
formed  hinge  is  "travelling"  or  not,  respectively.  However,  at  the  center  of 
the  bay  (x  ■  1)  the  slope  must  be  zero 


w^  (1,  t)  ■  0 


(A-56) 
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Due  to  symmetry,  the  shear  at  the  center  must  be  zero,  i.e., 

( 

(1,  t)  -  0  (A-57) 

Also  observe  that  the  non-dimensional  bending  moment  and  membrane  force 

n.  must  lie  on  the  yield  locus  when  the  structure  is  deforming. 

♦ 

From  beam  theory,  recall  that  the  bending  moment  distribution  of  a  clamped 

beam  has  a  maximum  negative  value  at  the  supports  (x  ■  0  and  x  ■  2) ,  while  it 

attains  a  maximum  positive  value  at  the  center  (x  ■  1).  Interpreting  this  in 

terms  of  the  square  yield  surface,  the  bending  moment  distribution  could  lie 

anywhere  along  the  DA,  AB,  BC  (including  corners)  portions  of  the  Tresca  locus. 

However,  we  also  know  that  due  to  the  compressive  nature  of  the  external 

pressure,  n^  must  be  negative.  This  limits  us  along  BA  and  the  portion  of 

AB  (Including  corners)  where  n.  Is  negative. 

9 

GENERAL  METHOD 


As  already  outlined  in  the  "PROBLEM  STATEMENT  AND  SOLUTION  METHOD”  section 
of  the  body  of  the  report,  the  approach  consists  of  solving  the  equilibrium 
equations  involving  two  stress  resultants,  m  and  n.,  subject  to  initial, 

X  $ 

boundary,  and  jump  conditions.  The  load  consists  of  a  non-dimensional  pressure 
load,  exceeding  the  collapse  load  pQ  (References  A-13,  A-14,  and  A-l)  which 
was  calculated  by  a  limit  analysis  and  a  rigid,  perfectly  plastic  material.  In 
the  process,  we  must  assume  (1)  a  kinematically  admissible  velocity  profile,  and 
(2)  the  portion  of  the  applicable  yield  locus. 


A-l  7 
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Based  on  these  assumptions,  we  obtain  ‘a  solution.  We  must  verify  its 
validity,  i.e.,  none  of  the  stress  resultants  must  violate  the  yield  surface. 
It  is  these  conditions  that  impose  certain  constraints  on  our  parameters,  and 
allow  certain  cases  to  emerge. 


RING 

STIFFENER 


RING 

STIFFENER 


NOTE:  GLOBAL  FRAME  OF  REFERENCE  (X.  Y.  Z)  IS  THROUGH 
LEFT  END,  WHERE  FIRST  RING  STIFFENER  IS  LOCATED. 


FIGURE  A-1.  CIRCULAR  CYLINDRICAL  SHELL  OF  LENGTH 
2L,  DIAMETER  2a,  SKIN  THICKNESS  h  =  2H 
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TABLE  A-l.  REGIONS  ON  SIMPLIFIED  TRESCA  SQUARE  (YIELD  LOCUS) 
Plastic  Stress  Resultants  Strain-Rate  Vector 


Region 

Components 

n4>  mx 

Inequalities 

Components 

w  w"/2c^ 

Equations 

Inequalities 

AD 

* 

-1 

_1  1  n*  <  1 

0 

X 

w  ™  0 

w”  >  0 

A 

-1 

-1 

* 

M 

X 

* 

w  >  0,  w"  >  l 

AB 

-1 

* 

-1  <  mx  <  1 

U 

0 

w"  -  0 

w  >  0 

B 

-1 

+1 

* 

U 

-v 

* 

w  >  0,  w”  <  i 
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CASE  A 


APPENDIX  B 


-  SHORT  SHELLS,  LOW  LOADING 


l  +  liPl1  +  l 

c  c 


(° 


<  c 


B-l 
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SHORT  SHELLS,  LOW  LOADING,  l+^x<p<l+^  WITH  c2  >  0 

c  <r 

In  this  section,  we  show  the  relevant  stress  resultants  m  and  n  , 

x  <t> 

displacement,  velocity,  and  acceleration  (w,  w,  w),  as  well  as  displacement 

at  rest  w(x,  t  ),  and  time  to  come  to  rest  t  in  Tables  B-l  and  B-2. 
o  o 

Figure  B-l  displays  the  kinematic  assumption  with  regard  to  the  velocity  field. 
The  analysis  proceeds  as  follows. 


The  non-dimensional  collapse  load  pQ  was  derived  by  Hodge  in  Reference  B-l. 
It  is  defined  as 


c 


We  must  solve  (Equation  A-ll  of  Appendix  A) 


— r  m"  +  n  +  p  -  w  =  0  (B-2) 

2c2  X  * 


subject  to  initial,  boundary,  and  jump  conditions  for 


p  >  p  >0  ( B— 3) 


and  p  *  0. 


B-2 
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Equation  (B-2)  must  be  solved  in  two  time  intervals. 

1.  For  0  £  t  £  1,  when  p  j4  0  and 

2.  For  1  _<  t  £  x  ,  when  p  ■  0, 

where  tq  represents  the  time  at  which  the  shell  comes  to  rest. 

SOLUTION  FOR  TIME  INTERVAL  0  <  t  <  1  ,  (p  J4  0) 

Referring  to  Figure  A-3,  assume  we  are  on  the  AB  side  of  the  Tresca  square. 
Using  Table  A-l,  this  translates  to  the  following  requirements  for  the  two  resul¬ 
tants  and  velocity 


n  -  -1  (B-4) 

-1  1  mx  1  1  (B-5) 

w"  -  0  (B-6) 

w  >  0  (B-7) 


B-3 
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Since  w"  -  0,  the  velocity  profile  can  only  be  linear  in  x.  Since  we  have  two 
initial  conditions  (zero  displacement  and  velocity  at  time  t  “0)  we  assume  a 
displacement  profile  of  the  form  (linear  in  distance  from  ring  stiffener) 


w(x,  t) 


[Ao  + 


*1 


T  + 


A2t2] 


(B-8) 


with  w  =  x[A^  +  2A2t] 


(B-9) 


w  **  2A2x 


( B-10 ) 


But 


w(x,o)  =*  0 


(B-ll) 


implies 


A  -  0 
o 


and  w(x,o)  *  0 


implies 


(B-12) 


(B-13) 


A 


1 


0 


(B-14) 


Therefore, 


w(x,  t) 


xA2t 


B-4 


(B-15) 
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I 

where  ^  f®  aot  a  function  of  time  or  distance  but  a  function  of  the  load 
level,  i.e., 


A2  "  A2(p) 


Also,  in  view  of  w  >  0 


a2  >  ° 


We 


now  replace  n.  *  -1  and  w  ”  xA-t  in  the  equilibrium  equation. 


and  solve  for  m 


m”  +  2c2(p-l)  -  4c2xA2  *  0 


m”  “  4c2A2x  -  2c2(p-l) 


Integrating  once 

m^  -  2c2A2x2  -  2c2(p-l)x  +  ^  (B-20) 

But  at  x  *  1,  m'(l,x)  ■  0,  i.e., 
x 

Bx  -  -  2c2A2  +  2c2(p-l)  (B-21) 

m'  -  2c2A.(x2-1)  +  2c2(p-l)(l-x)  (B-22) 
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Integrating  again, 


mx  *  +  2c2(p-D(x  -  iX2)  +  Dx 


3  “  c2(p-l)x2  +  2c2(p-l-A  )x  +  D 


At  the  end  support  (x  “  0),  m  ■  -1 


This  defines 


(B-23) 


D1  -  -1 


(B-24) 


Thus, 


mx  -  |c2A2x3  -  c2(p-l)x2  +  2c2(p-I-A  )x  -  1 


(B-25) 


At  the  middle  (x  -  1),  m  *  1. 

x 


This  determines  that 


A,  =*  -(p-1)  - 


2c' 


(B-26) 


B-6 
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an d,  finally,  the  bending  moment  profile  assumes  the  form 


"2  1  r  2 

2~(p  -1)  -  1  x3  -  c2  (p  -l)x2  +  j-(p-l)  +  3  x  - 


(B-27) 


and  the  displacement  w  becomes 


5?  Ix<p'1)  -  J1’2 


(B-28) 


while  the  velocity  w  is 


w 


2  XT 


(B-29) 


Now  we  must  require  that  |m^|  £  1  everywhere.  We  obtain  the  first  and 
second  derivatives  of  m^.  They  are 


8m  [2  1  9  9  (2 

~  -  3[j-(p-l)  -  ljx  -  2c  (p-l)x  +  j|-(p-l)  +  3 

-  6  ^(p-1)  -  lj  x  -  2c2(p-l) 


(B-30) 


(B-31) 


B-7 
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We  observe  that  one  of  the  roots  of  Equation  (B-30)  is  x±  -  1  (midpoint).  The 
other  root  is 


2 

c 


3  +  |  (p-1) 

,  i 

11  4 

r  2 

3 

n  t 

Lr(p-i)  - 1 

[  |r(p'1)  - 

(B-32) 


which  is  a  positive  quantity  since 

2 

p>p0*l  +  —  .  Asp  approaches  pQ  (from  the  right),  x  -*•<*>  . 

c  2 


This  meaus  that  if  we  choose  this  root  (x2)  to  exceed  or  equal  1,  in  the 
interval 


0  <  x  <  1 


Equation  (B-30)  will  always  be  non-negative,  i.e.. 


1 

3 


1  + 


4 


f^(p-l)  -  1 


>  1 


p<  i  +  -| 

c 


( B-33 ) 


B-8 
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We  need  to  examine  later  the  sign  of  the  second  derivative,  Equation  (B-31),  to 
determine  whether  the  yield  condition  will  be  violated.  We  started  the  analysis 
based  on  the  fact  that  for  plastic  deformations  to  take  place 


P  >  Pn  "  1  +  S 


means  that  ^-(p-1)  -  1  _>  0 


This  also  means  that  p  -  1  >  0 


At  the  left  end  (ring  support),  x  ■  0 


Y  m  ~  2c  (p-1)  <  0 


At  the  center  (x  ■  1) 


®  c  ^  2  2 

-  6  ~(p-l)-l  -  2c  (p-1)  -  c  (p-1)  -  6 

3xZ  <- 


Since  j-  (p-1)  -  1  >  0 


c  (p-1)  >  2 


(B-34) 


( B-35 ) 


(B-36) 


( B-37 ) 


(B-38) 


(B-39) 


B-9 


c2(p-l)  -  2  +  e 
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Set 


where  e  >  0 


Then 


(B-40) 


3 


-4  +  e 


x“l 


(B-41) 


At  the  two  end  points,  the  second  derivative  varies  from  -4  -  2e  (which  is 
always  negative),  to  -4  +  e.  The  function  will  have  a  maximum  if  its 
first  derivative  vanishes  and  its  second  derivative  is  negative  or  zero. 

We  must  have 


3 


6 


2 

2-(p-D-l 


x 


2c2(p-l)  £  0 


(B-42) 


V  /  _ 

c2(p-l) 

■  2 

3 

|-(P-1)  -  1 

for  all  0  <  x  <  1. 


At  x-1  -  4  +  e  <  0 


or  e  <  4 


B-10 


(B-43) 


(B-44) 


(B-45) 


i 
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( B-46 ) 


or 


c  (p-1)  “  2  +  e  <  6 


or 


p  <  1  +  ^ 


(B-47) 


2 

At  the  lowest  possible  load  p  ■  *  1  +  —  the  moment  resultant  assumes  the 

c 

form 


lower  ,  v  „  2  ,  . 
m^  (x,x)  *  -  2x  +  4x  -  1 


(B-48) 


while  at  the  upper  possible  load  (without  violating  the  yield  surface) 


P  "  1  +  \ 


(B-49) 


m^PPer  (X,T)  -  2x3  -  6x2  +  6x  -  1 


(B-50) 


Table  B-l  summarizes  the  obtained  results  so  far. 


SOLUTION  FOR  TIME  INTERVAL  l<T<T,(p-0) 

o 

At  time  t  *  1  the  pressure  ceases  acting  (p  *  0).  There  is  an  acting 
acceleration,  velocity,  and  displacement.  They  must  match  with  the  solution  in 
this  range.  Therefore,  at  t  ■  1 
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w  (x,l) 


|  (p-1)  -  ^  x  “ 
c  p“0 


')*?]*  +1] 


Furthermore,  the  moment  resultant  with  p  =  0  becomes 


"  2  "1  T  2' 

/  c  ,  .  3  1  22.  ,  c 

mx(x,l)  **  —  2  fix  +  cx  +  3  -  y~  x  - 


Assuming  the  range  AB  on  the  yield  surface, 


w  ’  «  0 


w  >  0 


to  follow  normality  requirements.  This  means  that  the  velocity  profile  must  be 
linear  in  x  (w”  =0). 


Assume 


w(x,t)  »  x[Bq  +  B^t  +  B2t  ] 


w  ■  x[B^  +  2B2t] 


w  =  2  B»  x 


(B-5 
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But  from  Equation  (B-51) 


(B-58) 


•  h  [r  +  ^ 


(B-59) 


Therefore, 


w  -  x  Bi  -  +  l]r 

c  '  ' 


(B-60) 


Again,  using  normality  (w  >  0) 


’liitr*1] 


(B-61) 


Furthermore, 


p  > 1  +  r 


(B-62) 


and,  hence, 


B1  >  f  p 


(B-63) 


B-13 
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Taking  the  lowest  value  of  for  which  normality  is  satisfied, 


w 


( B-64 ) 


Therefore, 


w  >  0 


(B-65) 


gives 


3 

2 


It  >  0 


(B-66) 


or 


(B-67) 


When  time  t  becomes  tq,  velocity  is  zero  if  the  shell  comes  to  rest. 


Consequently,  t 


£_ 


Also, 


w 


Bo  +  2pT 


(B-68) 


B-14 
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Since  displacements  must  match  for  time  t  ■  1 


2 

^-(p-1) 


yielding 


(B-69) 


B  » 
o 


and 


w  =  x 


+ 


(B-70) 


(B-71) 


The  displacement  at  time  t  ■  tq  can  be  calculated  by  substituting  the 

value  of  r  .  Hence, 
o 


We  must  solve  again 


m 


+  % 


+  P 


w 


0 


(B-73) 


with  p  *  0,  n. 


-1  and  w 


2 

£_  + 
2  L2 


(B-74) 


x  + 


1 


(B-75) 


B-15 
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m" 

z 


(B-76) 


Integrating 


m* 

x 


+  2c  x  +  E, 


(B-77) 


However,  this  must  match  the  solution  of  time  t  *  1,  i.e.. 


m^(x,l)  -  -  3^  +  l]x2  +  2c2x  +  ^  - 


-  3 


+  1 


a2+  2c2x  +  [s  -  §-] 


(B-78) 


Hence, 


E 


1 


3 


(B-79) 


Integrating  again 


m 


x 


2 

c 

T 


+  1  x"  + 


2  2 
C  X 


x  +  E 


2 


(B-80) 


This  must  match  the  solution  at  time  r  ■  1,  i.e.,  m  ~  1 


B-16 


NSWC  TR  86-328 


Therefore, 


m 

x 


c2  1  1 

-  F  *  y* 


2  2 
c  x 


1 


We  observe  again  that 


(B-81) 


mx  (0,t)  -  -1  (B-82) 

m  (1,t)  -  1  (B-83) 

x 

Again  the  yield  condition  must  not  be  violated  along  any  portion,  which  means 
that  -1  <_  _<  1.  The  moment  resultant  is  an  increasing  function  from  x  ■  0  to 

x  ■  1  and,  for  it  to  stay  that  way  without  exceeding  the  maximum  absolute  value 
of  1,  the  first  derivative  must  be  non-negative  in  the  x  range  (0,1).  Further¬ 
more,  near  x  -  1,  the  second  derivative  must  be  non-positive  for  an  increasing 
function. 


The  first  condition  is 


(B-84) 


Since 


1 


>  0 


(B-85) 


B-17 


NSWC  TR  86-328 


the  quantity  in  curly  brackets  must  be  non-positive.  For  this  to  be  so,  the 
value  of  x  must  lie  within  the  range  of  the  two  roots  of  the  quadratic.  Since 
by  inspection  one  of  them  is  ■  1,  the  other  must  necessarily  be 


x 


2 


(B-86) 


Now  x  will  always  be  in  the  non-negative  range  0  £  x  £  1,  which  means  that  if  we 
make  x^  negative  or  zero,  x  will  always  satisfy 


Therefore, 


<  0  (B-87) 

gives  c2  £  6  (B-88) 

The  second  condition  is  always  satisfied  for  points  x  *  1  -  6  (when  1  >  6  ^  0, 
and  6  is  a  small  quantity) 


-  -  6  +  lj(l-fi)  +  2c2  _<  0 

x"l-6 


(B-89) 


B-18 
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or  when 

(B-90) 

2 

For  small  c  (which  is  the  current  case),  6  must  be  less  than  1  (as  assumed). 

Finally,  Figure  B-l  summarizes  the  velocity  profiles  for  this  case,  and 
Table  B-2  summarizes  the  used  quantities. 
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TABLE  B-1.  SUMMARY,  SHORT  SHELLS,  LOW  LOADING,  1 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE  j 

SHORT  (O  <  C2<6) 

CONDITIONS 

0  <  T  <  1 

MOMENT 

RESULTANT 

mx  (x  ,r )  *  £-£-  (p  -  1)  -  ijx  3  -  C2  (p  -  1)x  2  +  £_2?(p-  1)  +  aj  x  -  1 

MEMBRANE 

RESULTANT 

=  -1 

DISPLACEMENT 

VELOCITY 

W  <x. T )  «  -^2 [c2(p-11- 2> 

ACCELERATION 

w(x,r)  » -^2|c2<p- 1) -2jx 

TIME  , 
ro 

N/A 

DISPLACEMENT 
AT  REST 

HAS  NOT  COME  TO  REST  YET 

B-2I 
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TABLE  B-2.  SUMMARY,  SHORT  SHELLS,  LOW  LOADING  ,  2 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

short  |o  <  c2<  ej 

CONDITIONS 

1  <  r 

Ho 

MOMENT 

RESULTANT 

mx  (x  ,r  )  =  - 

f4-' 

MEMBRANE 

RESULTANT 

1 

II 

s- 

c 

DISPLACEMENT 

w  (x,  r)  -  £ 

_ 

VELOCITY 

ACCELERATION 

W(x,r)  =  --^(^+1)x 

TIME  . 

o 

_  p  c2  „ 

T°  Po  (c2  +  2) 

DISPLACEMENT 
AT  REST 

B-22 


NSWC  TR  86-328 
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CASE  B 


APPENDIX  C 


-  LONG  SHELLS,  LOW  LOADING 


2 

— j  <_  p  <  1  + 

c 


«) 


C-l 
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LONG  SHELLS,  LOW  LOADING, 


\  <  P  <  1  + 
c 


WITH  c2  >  0 


Figure  C-l  displays  the  assumed  kinematics  for  the  velocity  distribution  for 
times  0  £  t  £  1  and  1  £  t  £  tq.  Tables  C-l  through  C-5  summarize  the  results  we 
are  about  to  obtain.  Observe  that  there  are  two  intervals  of  interest,  depending 
on  whether  o  £  t  £  1  and  1  £  t  £  xQ.  Furthermore,  the  last  interval  is  subdi¬ 
vided  depending  on  whether  the  point  of  interest  x  is  to  the  left  (o  £  x  £  u)  or 
to  the  right  (u  <  x  <  1)  of  the  travelling  hinge. 


When  the  pressure  loading  is  restricted  in  the  range 


1+ 


^  <  p  <  1  +  ^ 

c  c 


(C-l) 


the  moment  resultant  m^  of  the  previous  analysis  [Equation  (B— 27 ) ]  for  o  £  x  £  1 

does  not  violate  the  yield  surface  requirements  -1  £  m^  £  1.  In  the  second 

2 

range,  however  (1  £  t  £  tq)  [Equation  (B— 81) ] »  when  c  >  6,  we  obtain,  setting 

2 

c  ■  6  +  c  and  e  >  0  and  e  being  small 


m^ ( x , r ) 


4  +  |  x3 4  +  (6+e)x2  -  |  x  -  1 


(  C— 2  ) 


It  is  quite  obvious  that  if  x  ■  ^  (i.e.,  near  a  point  where  x  ♦  +  0),  the 

yield  condition  is  violated  since 

C-2 
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(3 


T) 


[4e3+90e2+864e+3456]  <  _  ± 
54[8+e]2  " 


(C-3) 


Table  C-l  summarizes  the  result  for  0  <  t  <  1.  We  just  observed  that  when 
2 

c  >6,  the  yield  condition  is  violated  near  x  ■  +0. 

Consider  the  possibility  of  two  regions  on  the  yield  surface: 


1.  Region  AD  for  point  x  such  that  0  <_  x  £  u 

2.  Region  AB  for  point  x  such  that  u  <  x  £  1 

where  u  is  the  point  along  the  shell's  length  where  the  regions  change.  It 
actually  represents  the  position  where  a  hinge  circle  develops. 

REGION  AD 


The  hinge  in  question  is  travelling  because  it  changes  position  with  time. 
Its  initial  location  must  be  at  the  ring  support  at  time  t  ■  1.  Along  AD,  the 
strain  rate  vector  must  satisfy 

1.  w  •0  (no  velocity)  (C-4) 

2.  w"  -  0  (C-5) 

3.  »x  -  -1  (C-6) 


C-3 
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and  equilibrium 


— m"  +  n.+p-w“0 
2c2  x  ♦ 


(C-7) 


(See  Figure  A-3.) 


Since  w  *  w  *  0  and  at  time  t  ■  1  the  pressure  has  ceased  acting  (p  *  0) 


n  -  0  (C-8) 

v 

and 

w  -  cx  (C-9) 

The  deformation  is  rigid  plastic  for  0  <  x  <  u.  This  means  that  all  points 
0  <  x  <  u  have  no  motion  for  1  <_  x  <_  tq,  i.e.,  they  do  not  deform 
further  than  the  deformation  they  acquired  at  time  x  ■  1,  which  is 


w(x,l) 


2 

— (p-D-l 


2 


x 


(C-10) 


REGION  AB  FOR  u  <  x  <  1 


The  conditions  the  strain  rate  must  satisfy  are 


w" 

-  0 

(C-ll) 

w 

>  0 

(C-12) 

% 

-  -1 

(C-13) 

C-4 
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where 


-1  <  a  <1 
—  *  — 


while  equilibrium  gives: 


— -  m"+n  +  p  -  w  ■  0 
2c2  x  * 


(C-14) 


Since  the  second  derivative  of  the  velocity  with  respect  to  position  must  vanish 
it  must  be,  at  the  most,  linear  in  distance.  Assume,  therefore, 


and 


w  -  A  x  +  B 


w  «  A  x  +  B 


(C-15) 


(C-16) 


Replacing  Equation  (C-16)  in  the  equilibrium  Equation  (C-14)  and  setting  p  ■  0, 
we  have 


■4  m"  -  1  -  (A  i  +  S)  •  0 

2c  X 


(C-17) 


Integrating  once 


m^  ■  2c2x  +  2c2  ^|a  x2  +  §x^  +  c^ 


(C-18) 


C-5 
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Integrating  again 

22  2  flv  3  .  u  2"!  , 

mx  -  e  x  +  c  [3A  x  +  fi  x  J  +  C]x  +  c2 


(C-19) 


At  the  two  ends  x  ■  u  and  x  ■  1  the  moment  resultant  becomes  -1  and  1, 
respectively. 


mx  -  c2x2+  c2  x3+  fi  x2]  +  c^x  +  c2  (C-20) 

At  x  ■  u  mx  *  -1.  Therefore, 

c2  £^Au3  +  (l+fi)u2J  +  c^u  +  c2  -  -1  (C-21) 

Continuity  of  shearing  force  (m^  *  oj  at  x  ■  u  yields 

m^(u)  -  c2[a  u2  +  2(1  +  fi)  u]  +  cx  -  0  (C-22) 


At  the  middle  (x  ■  1)  we  must  have 

.,(1)  -  i 

»'(1)  -  0 
X 

i.e. 

c2[A  +  2(1  +  B)]  +  ci  -  0 

c2  +  (1  +  fi)J  +  ci  +  c2  ■  1 

C-6 


(C-23) 


(C-24) 


(C-25) 

(C-26) 


1 


a 


1 
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cx  -  -  c2[X  +  2(1  +  8) ] 


c2  -  !  +  e  +  (1  + 


Substituting  Equation  (C-27)  for  in  Equation  (C-22)  we  get 


-2(1  +  8)  -  -2(1  +  8) 

(«  -l) 


n+sr 


for  u  +  1 


Substituting  and  in  Equations  (C-21)  and  (C-22)  we  get 


sl<»rl.>(;2  +  .u-2)x  +  c2(1  +  6)(u2  .  2  u  +  1)  ■  -! 


However,  from  Equation  (C-29) 


1  +  8  -  -  j(u  +  1)X 


and  after  algebra  we  obtain 


A  -  -  ^  _L 

c2  (1-u)3 


1+8.6  (u±l) 

c  (l-u)J 


(C-27) 


(C-28) 


(C-29) 


(C-30) 


(C-31) 


(C-32) 


(C-33) 


C-7 
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This,  in  turn,  defines  and  from  Equations  (C-27)  and  (C-28).  They  are 


(1-u)' 


(C-34) 


„  _  ,  _  2(l-3u) 
2  (1-u)3 


(C-35) 


The  two  solutions  [Equation  (C-15)]  in  the  intervals  0  <  x  <  u  and  u  <  x  <  1 
must  have  continuous  velocities  at  x  -  u,  i.e., 


w  =  A(t)  u  +  B(t)  ■  0 


(C-36) 


Differentiating  Equation  (C-36)  again  with  respect  to  time  we  obtain 


Au  +  Au  +  B“0 


(C-37) 


i  .  _  ,v  .  _  6  c4  1 

A  u  -  -(A  u  +  B)  *  — r-  i - - 

c  .  (1-u) 


(C-38) 


A  -  1  -  6  1 

u  2  ...  s 2 
c  u(l-u) 


(C-39) 


Differentiating  Equation  (C-39)  with  respect  to  time  we  get 


V  .  U_  6_  -(u(l-u)2)1 

'u2  c2  [»Cl-»)2]2 


ii  6  u  12  1 

“.22  .2,.  x2  “  2  *3 

u  c  u  (1-u)  c  (1-u) 


(C-40) 


C-8 
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He  already  obtained  [Equation  (C— 32) ]  that 


X-.12 

cr  (1-u)' 


Equating  the  two,  we  obtain  the  following  differential  equation 


Since 


m-S-l-o 

u  Lc  (1-u)2  J 


*-Sx  f  --i-,  il*0 

u  (1-u) 


we  must  necessarily  have 


H-o 

u 


ii  ■  0 


u  +  0 


This  means  that  the  location  of  the  hinge  circle  Is  linear  In  time,  x 


U  *  Ej_  T  +  E2 


u  ■  Ei 


ii  -  0 


(C-41) 


(C-42) 


(C-43) 


(C-44) 


(C-45) 


(C-46) 


(C-47) 


(C-48) 


(C-49) 


C-9 
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Velocities  must  match  at  time  t  ■  1.  [See  Equations  (C-36)  and  (B-29).] 
Therefore, 


-(p  -  1)  -  1  tx  -  A(l)x  +  B(l) 


(C-50) 


i.e.,  we  get  two  initial  conditions  for  A,  B 


[r(p"1) ' 


(C-51) 


0  **  B(l) 


(C-52) 


From  Equation  (C-36) 


B  ■  -  Au 


(C-53) 


and  from  Equation  (C-43) 


s  _  6_  c~_  1  1 

2  6  "  .2  u 

c  L  (l-u)  J 


(C-54) 


Therefore,  by  Equation  (C-50) 


1 6 c_ 

2  6 


l-u)  J  ,-i  c 


(C-55) 


C-10 
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with 


u(t)  -  EjT  +  E2 
u(t)  -  Ej^ 
and  u(l)  -  0 


Equations  (C-58)  and  (C-56)  imply  E^  ■  -  E^,  i.e., 
u(t)  -  E-^r-l) 


Replacing  Equations  (C-59)  and  (C-60)  at  time  t 

nht-a<  n 


(C-56 

(C-57 

(C-58 

(C-59 

(C-60 

1  in  Equation  (C-55)  we 

(C-6i: 

(C-62! 

(C-63! 


C-ll 
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Therefore,  by  Equation  (C-53)  the  velocity  field  for  u  <  x  <  1  is 


w  *  A(t)x  +  S(t)  ■  A(t)(x  -  u) 


However,  using  Equations  (C-54),  (C-62),  and  (C-63) 


x  _  6  1 

z 

c 


(1-uV 


(x  -  u) 


where  u  and  u  are  given  by  Equations  (C-62)  and  (C-63),  respectively. 


Replacing  the  values  of  c^,  c^,  X,  and  8  in  Equation  (C-20)  we  obtain 


mx(x,t)  -  c 


i.  _J: _  r3  +  6.  (1-Hi)  r  2 

c  (1-u)-3  c  (1-u  V 


12  .  ,  2(l-3u) 

- ~  U  X  +  1 - - - 7T- 


(1-u)' 


(1-u)' 


1  +  __i — L  4X3  +  6(l+u)x2  -  12ux  -  2(l-3u)l 
(1-u)-3  L  J 


Also 


%  “-1 


1  <  T  <  T 
—  —  O 


(C-64) 

(C-65) 


(C-66) 

(C-67) 

(C-68) 


C-12 
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Equation  (C-66)  yields 


0^(1, t)  -  1 


®  (u,t)  -1+2 


w  >  0 


,(u-l)J  _ 
■ - 5 


w"  -  0 


Equation  (C-71)  Implies 


6  1  c 


(C-69) 


(C-70) 


(C-71) 


(C-72) 


c2  u  |_6  (1-u)2 


(x-u)  >  0 


(C-73) 


with  u  given  by  Equation  (C-63)  and,  hence,  u  given  by  Equation  (C-62) 


3  I  §-(p-l)-l| 


(C-74) 


Since  u  >  0  (c2  6) ,  Equation  (C-73)  can  be  rewritten  as 


5,  ■  (f  ^  --  k2.a=“4^1  <*-»>  i  o 

c  3  U-(p-l>-ll  >■  *a_u)  J 


(C-75) 


C-13 
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This  implies  that 


2 

c 


(1-u)2 


6  >  0 


(C-76) 


or 


1 1-u  | 


(C-77) 


Since  1  j>  u  >  0,  Equation  (B-77)  reduces  to 


u  >  1 


c 


(C-78) 


At  time  x  ■  tq  (>^  1)  the  motion  stops 


w(x,tq)  *  0 


(C-79) 


i  •  e  m  f 


*(To)[X-“(ro)] 


-  0 


(C-80) 


or  since  in  general  x  j*  u 


A 


-  0 


(C-81) 


or  by  Equation  (C-54) 


(C-82) 


C-14 
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Since  in  general,  u(t0)  ■ 


(«*-•) 


>  0 


by  Equation  (C-82) 


yielding 


i 


(C-83) 


i 


u 


(C-84) 


Replacing  t  *  tq  in  Equation  (C-63)  and  equating  it  to  Equation  (C-84),  we  solve 

for  t  from 
o 


/S'  _  1 

(c2-*) 

c  3 

■2  1 

§-(p-l)-lJ 

1 

[2 

3 

If-Cp-iM 

-  1  +  _J 

c  (c+/S) 

or  if  we  set 


(C-85) 


(C-86) 


a 


3 


§-(p-l)-l 


(C-87) 


C-15 
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(C-88) 


and  Equation  (C-65)  for  the  velocity  field  can  be  rewritten  as 


w  ■ 


1 

(1-u)2 


(x-u) 


1 

rTT+oT-arJ 


+  (»"1) 

[ (l+a)-at] 2 


Also  observe  that 


and 


and 


(C-89) 


(C-90) 


(C-91) 


(C-92) 


Integrating  Equation  (C-89)  and  applying  the  boundary  condition  that  [Equation 
(B-28)  ] 


w(x,l) 


(C-93) 


C-16 


NSWC  TR  86-328 


we  have 


w<x»t)  ■  —j  £|-(p-l)-ljx  +  j * w 


(C-94) 


Observing  that  (see  Tables  C-2  and  C-3) 


—  -  -log  |  (l+a)-<xr  I  +  "1-  - — 

dt  a  e  „  n*«- 


1 _  +  (x-1) 


a  [1+a-ar] j  [1+a-ax]  [1+a-at]2 


(C-95) 


Equation  (C-94)  becomes 


■  £1  [&Pl>-l]«  *  ^  (I)  _/|^t«-a(x-l)l  - 

jdT‘^  [v^-1’-1]1 


S®  jr[(x+a)(T-1)  -|(t2-i)]  + 


(C-96) 


C-17 
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or  using  Equations  (C-90),  (C-91),  and  (C-92)  the  displacement  at  time  tq 
[i.e.,  when  w(x,x  )  ■  0]  becomes: 


(C-97) 


Observe  that  x  *  u  only  at  time  x  «  1.  At  that  time,  x  ■  u  ■  0  and  Equations 
(C-96)  and  (C-10)  give  the  same  answer. 


Equation  (C-96)  can  also  be  written  in  terms  of  u  as  follows: 


w(x,x) 


(x+q)u 

a 


(x+1) 


+ 


1 

a 


logell-u| 


+ 


(l-x)u) 
(1-u)  1 


(C-98) 


The  acceleration  w  is  given  by  replacing  Equations  (C-32)  and  (C-33)  in  (C-16) 


Ax  +  B 


K  — - — »[  (l+u)-2x]  -  1 
c  (i-u r 


(C-99) 


C-18 
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If  we  further  substitute  u  from  Equation  (C-63)  we  get 

| 

- i - ,[1+o(t-1>-2x]  -  1  (C-100) 

c  [1+a-aTj 

I 

Figure  C-l  summarizes  the  velocity  profiles  employed  for  long  shells  (c2  >  6) 
and  low  loading  |l  +  ^j<p<l  + 

I 

Tables  C-4  and  C-5  are  a  summary  of  all  the  derived  quantities  for  the  same  case. 


t 

i 

I 

i 


C-19 
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TABLE  C-1.  SUMMARY.  LONG  SHELLS,  LOW  LOADING,  1 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

LONG  (c2  >  s) 

CONDITIONS 

0  <  T  <  1 

MOMENT 

RESULTANT 

mx  (x.r)  x3  —  C2  (p  —  1 )  x2  +  |^(p-1)  +  3jx-1 

MEMBRANE 

RESULTANT 

v  -1 

DISPLACEMENT 

w(x,r)  «=  _2_  \Sl.  (p  —  1)  —  llx  t2 

VELOCITY 

w  (x,t)  ■  jtjjc^tp-U-sJw 

ACCELERATION 

w  (x,r)  -  JL  C2(p-1)-2|x 

2c2L  J 

TIME,. 

To 

N/A 

DISPLACEMENT 
AT  REST 

HAS  NOT  COME  TO  REST  YET 

C-21 
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x  x 


TABLE  C-4.  SUMMARY,  LONG  SHELLS,  LOW  LOADING,  2 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

LONG  jc2  >  6j 

CONDITIONS 

0  <x  <u  u  —  Od(x  — 11  Ot  ~  (c26)  ,1<t<t 

u  •  01  ,c 2  . 

3[|-(p-1)-l]  C 

MOMENT 

RESULTANT 

mx  -  -1 

MEMBRANE 

RESULTANT 

11 

o 

DISPLACEMENT 

w  =  JL  f£?(p-i)-il  x 

2C2  l  2  ' 

VELOCITY 

w  =  0 

ACCELERATION 

vy  =  0 

TIME  T 

o 

O 

II 

DISPLACEMENT 

AT  REST 

w(x,tq)  =  _L_^  f  £f  (p-l)-lj  x  FOR  0  **  *  <  a  (r„  —  1) 

2C2  1  2  ' 

C-24 
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TABLE  C-5.  SUMMARY,  LONG  SHELLS,  LOW  LOADING,  3 


.TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

LONG  ^C2  >  6^ 

LOW  (\  +  _L<p  <  1 

V  c2  C2' 

CONDITIONS 

(C2  -  6) 

O'  “  .  .  1  <  r  <  u  <  r  <  1 

3[^(p-1)-l] 

c 

MOMENT 

RESULTANT 

mx  (x,t)  »  1  +  1  4x3  +  6(l +u)»2  -12ux  -2(1  -3u)l 

(1  -  u)\  J 

where  u  *  ot(r  —  1) 

MEMBRANE 

RESULTANT 

n  .a _ 1 

'p  1 

DISPLACEMENT 

-a(T  -1)1  +  ^  _x^.  11 

a  (1  +a-ar  )Jj 

VELOCITY 

w(x,r)  -  JL(1)  [.£?-  — \ —  1  (x-u) 

c2  W  l  6  (i  _  u)2  J 

ACCELERATION 

W  (x,  T  )  = 

/  6  \  [1  +a(r—  1)  —  2  x  ] 

\C2/  [l+a-ar]3 

TIME  T 

o 

T  «  1  +  l 

—  (p  -  1)  -  l] 

2  J  »  i  +  1  U  _Vf\ 

c(c+v<r)  a  \  c  f 

DISPLACEMENT 
AT  REST 

W(X,T0)=  2  [— (p-  1)  -  ll  X  + 

2C2  '  2  ‘ 

+l/± 

a2  \c2 

C-25/C-26 
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SHORT  SHELLS,  HIGH  LOADING , ^p  >  1+  ^  ,  0  <  c2  <  6 


Figure  D-l  gives  the  assumed  velocity  profile,  which  will  be  justified  in 
this  section. 


Tables  D-l,  D-2,  D-4,  D-5,  D-6,  D-7,  and  D-8  summarize  the  solution  in  all 
intervals.  Tables  D-4  and  D-5  summarize  the  solution  in  the  interval  1  £  t  _<  r', 
while  Tables  D-6,  D-7,  and  D-8  for  r *  <  t  <  t  . 

First,  we  observe  that  in  view  of  Equation  (B-42)  or  (B-47),  the  solution 
for  short  shells  (portion  AB  on  the  yield  surface)  in  the  range  0  £  t  _<  1  cannot 
be  used,  as  it  predicts  a  moment  resultant  mx  outside  the  (-1,1)  range  of  the 
Tresca  square. 

From  the  point  of  view  of  time  intervals,  it  turns  out  that  we  must  consider 
three  intervals: 

1.  STAGE  1  (0  _<  t  £  1)  for  times  during  which  the  excess  pressure  p  is 
acting. 

2.  STAGE  2  (1  <_  x  £  t’)  for  times  during  which  the  pressure  load  has  been 
removed  (p  ■  0),  but  motion  continues  (w(t')  j*  0),  t’  will  be  defined  there. 


D-2 
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3.  STAGE  3  (t'  £  t  _<  t  )  for  times  during  which  the  pressure  p  has  been 
removed  (p  ■  0)  and  the  shell  comes  to  rest  (w(tq)  *  0). 

STAGE  1 


We  consider  two  regions  as  a  possible  assumed  profile  on  the  yield  surface: 

(1)  for  points  such  that  0  <  x  <  uq  on  AB  (Figure  A-3)  and  (2)  for  points  such 

that  u  <  x  <  1  at  corner  A  (Figure  A-3).  Other  combinations  were  considered 
o  — 

and  eliminated  by  Hodge. ®  ^  They  violate  the  -1  _<  m  <_  1  requirement. 

1.  First  interval  0  <  x  <  ufl  (AB).  The  requirements  on  the  strain  rate 
vector  are  such  that 


w"  -  0 

(D-l) 

w  _>  0 

(D-2) 

“-1 

(D-3) 

Equations  (D-l)  and  (D-2)  together  yield  that  w  is  linear  in  distance  x,  such 
that 


w  *  Ax 


(D-4) 


with  A(t)  >  0 


(D-5) 


Therefore, 


w  ■  A(t)x 


(D-6) 


D-3 
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Equilibrium  requires 


m 


+  n. 


+  p  -  w 


0 


(D-7) 


or  using  Equations  (D-3)  and  (D-6) 


m” 

x 


+  2c 


(-1+p)  -  2c 


0 


CD-8) 


Integrating  twice  Equation  (D-8) 


m' 


x 


2c 


(D-9) 


m 


1  2 
3  C 


3  2  2 

Ax  -  c  (p-l)x  +  C..X  +  D, 


(D-10) 


At  x  ”  0,  mxC°)  =  "1>  yields 


D 


1 


*  -1 


(D-ll) 


Hence, 


m 


1  r2 

I  C 


A  x~ 


-  c2(p-l)x2  +  C. 


x  -  1 


(D-12) 


To  proceed  further  we  need  to  consider  the  second  region  uq  <  x  <_  1  and 

then  examine  both  solutions  at  x  ■  u  . 

o 
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2.  Second  interval  uq  <  x  £  1  (point  B  on  yield  surface), 
requirements  are 


w  >  0 


w"  <  0 


The  strain  rate 

(D-13) 

(D-14) 


and 


n.  **  -1 


(D-15) 


m  ”  1 
x 


(D-16) 


together  with 


m”+n  +p-w-0 
2c2  X  * 


Substituting  Equations  (D-15)  and  (D-16)  in  Equation  (D-17) 


w”n+p  =  p-  l 
$ 


and  integrating  once 


w  ■  (p-l)r  +  D.. 


D-5 


(C-17 ) 


(D-18) 


(D-19) 
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Since,  however,  w(x,o)  *  0,  “  0. 

Therefore, 


w  -  (p-1)t 


Integrating  again 


(D-20) 


w  *  jCp-1)t2  +  (D-21) 

where  E.  will  be  determined  by  matching  displacements  of  x  -  u  from  the  first 
A  o 

interval  (0  <  x  <  u  )♦ 

—  o 

Since  in  the  second  interval  (u  <  x  £  1)  m  -  1,  it  must  match  m  of  the 

^  x 

first  interval  (0  <.  x  <  uq)  ,  i.e,  by  Equation  (D-12) 

I  c2  A  -  c2(p-l)  u2  +  C1uq  -  1  -  1  (D-22) 

Furthermore,  the  shearing  forces  must  match,  i.e, 

-  [c2  Ax2  -  2c2(p-l)x  +  CJ  * 

X  ■  u 

0 

c2  a  u2  -  2c2(p-l)uQ  +  C1  -  0  (D-23) 


D-6 
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Since  the  hinge  may  be  stationary,  there  may  be  a  discontinuity  in  slope. 
However,  displacements,  velocities,  and  accelerations  must  be  equal  at  x  ■  uq, 
i.e.,  from  Equations  (D-6)  and  (D-18) 


A(t)uq  “  p-1 


(D-24) 


Therefore, 


A(t )  *  ^ 


(D-25) 


and  by  Equations  (D-25)  and  (D-6) 


w(T)  -  x 

o 


(D-26) 


represents  the  acceleration  in  the  first  interval  (0  £  x  <  uq)  .  By 
Equations  (D-4)  and  (D-19) 


A  uq  **  (p-l)x 


(D-27 ) 


i.e. , 


A(t )  “  x 

u 

o 


(D-28) 


D-7 
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This  means  that  u  is  a  constant  and  in  0  <  x  <  u 
o  —  —  o 


w(x,t)  ■  Ax  ■  -  t  x 

o 


and 


w(x,t) 


1  (p-1) 

2  u 

o 


2 

T  X 


(D-29) 


(D-30) 


because  w(x,o)  ■  0. 


At  this  point  we  consider  that  at  x  *  uq  the  displacements  must  agree.  By 
Equations  (D-30)  and  (D-21),  *  0.  Hence,  in  the  second  interval  (uq  <  x  _<  1) 


w  “  i(p-l)x2 


CD-31) 


Furthermore,  replacing  Equation  (D-25)  in  Equation  (D-12)  we  get 


for  (0  <  x  <  u  ) 

—  —  o 

Substitute  A  from  Equation  (D-25)  in  Equation  (d-23)  to  obtain 
C1  “  c2  (p-1)  uq 

D-8 


(D-32) 


(D-33) 
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Substitute  the  values  of  from  Equation  (D-33)  and  A  from  Equation  (D-25)  in 
Equation  (D-22)  to  determine  u 


o  2,  .  N 

c  (p-1) 


In  Equation  (D-32)  the  moment  resultant  can  be  written  in  an  alternate  form 


m^  (x,t)  -  2 


k '  ')' 


Therefore,  Tables  D-l  and  D-2  display  the  solutions  for  times  0  <  t  <  1  in 
the  two  intervals,  o  _<  x  <_  uq  and  uq  <  x  £  1.  Figure  D-l  gives  information  on 
the  velocity  profiles  for  times  o  <  t  <  1  as  well  as  1  £  t  £  t',  which  will  be 
analyzed  next. 


Note  once  more  that  since 


.  S_  ('*-  -  i\2  > 

,x  uo  V*.  J  - 


3  "x  .  12  /« 

,  2  2  (u 

3x  u  \  o 


-  1  <  0 


The  function  m^  is  uniformly  increasing  in  the  interval  0  £  x  £  uq.  This 
is  to  be  considered  later  (long  shells,  high  loading). 
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STAGE  2 


Pressure  load  p  is  removed,  p  ■  0,  where  t'  is  a  time  that  will  be  defined 
later.  We  must  solve  the  equilibrium  equation 

-TP  m'+p+n.  -w“0  (D-38) 

2c2  X  * 

with  p  ■  0,  with  initial  conditions  as  to  velocity  and  displacement  given  by 
Stage  1.  We  must  also  satisfy  the  conditions  imposed  by  the  normality  require¬ 
ments  of  the  strain  rate  vector  to  the  yield  surface.  We  assume  that  the  same 
two  ranges  on  the  yield  surface  apply,  i.e., 

1.  First  interval.  For  0  <  x  <  u,  we  are  along  AB. 


n.  -  -1  ( D-39 ) 

9 

-1  <  m  <1  (D-40) 

—  x  — 

w  >_  0  (D-41) 

w"  -  0  (D-42) 


and 


w(l) 


(?zi2  x 

u_ 


(D-43) 


D-10 
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u2-^- 

c2(p-l) 


CD-44) 


i 


w(l) 


_  (p-1) 

u 

o 


(D-45) 


w(x,l) 


1  <£li>  x 

2  u  * 

o 


(D-46) 


mx(x,l) 


&  - 


+  1 


(D-47 ) 


2.  Second  Interval. 
Figure  A-3). 


m  -  -1 
x 


w  >  0 


w"  <  0 


For  u  <  x  <  1  (at  corner  B  of  yield  surface;  see 


(D-48) 

(D-49) 


(D-50) 


(D-51) 


and 


w(l)  -p-1 


(D-52) 


u  2  -  <£li>  (D-53) 

o  u 

o 


D-ll 
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w(l)  -  p  -  1 

w(x,l)  “  ~(p-l) 

1.  First  Interval  0  £  x  £  u 
(1  <  t  <  t’) 


Take 


w(t)  -  (p-x)^ 

Then  by  differentiating  with  time  t 


Solving 


(D-54) 


(D-55) 


(D-56) 


(D-57) 


(D-58) 


(D-58) 


(D-59) 


D-12 


But 


Hence , 


Also  at  x 


(p— r)ux  - 


+  D 


1 


ix(o)  *  -1  gives  ■  -1 


c2x2  -  c2ux  +  c2(p—r)ux  - 


(i ♦ 


-  i 


■  u,  m  •  1  yields  the  following  ordinary  differential  equation 


2  2,  v  .  _  ,  1  2  2 

•=•  C  (p— r)uu  ■  2  +  -r  C  U 


Ifi 
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or 


2  (p-x)udu  ,  2_  ,  1  2 

3  dr  2  +  3  U 

c 


(D-66) 


dr  2udu 


(D-67) 


Hence  integrating  we  get 


-  logelp—r| 


+  C, 


(D-68) 


But  at  t  “  1 


2 

u  o - 

°  C  (p-1) 


(D-69) 


Therefore, 


C 


1 


(D-70) 


and,  finally, 


2 


u 


6t 

c2Cp~t) 


(D-71) 


D-14 
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Using  the  differential  Equation  (D-65)  in  Equation  (D-57),  the  acceleration  can 
be  rewritten  as 


(D-72] 


However,  using  Equation  (D-71)  we  get 


w 


(D-73) 


(D-74' 


which  makes  acceleration  linear  with  distance  for  points  0  £  x  £  u.  We  also 
note  that  it  is  a  "deceleration"  (negative  sign)  since  the  load  was  removed 
while  in  Stage  1  (0  £  t  £  1)  [Equations  (D-6)  and  (D-25)]. 

w(x)  “  (D-75] 

uo 


represents  a  positive  acceleration. 


Using  Equation  (D-65)  in  the  moment  resultant  Equation  (D-64),  we  obtain 


mx(x,r)  -  -  i(2  +  c2u2)(i)3  +  c2u2(i)2  +  (3  -  ^)(f)  -  1  (D-7 


which  becomes  1  at  x  ■  u. 
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Therefore,  for  0  £  *  £  uQ 


,  3/2  c2\  2  ,2  (3  c2u| 

■i’-Trs  +  i-J1  +2cx  +  ^-t j 


and 


2 


One  of  the  roots  of  *  o  is  "  u.  Thus  the  other  one,  p 
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Siace  m^(x)  ^  0  and  u,c  >  0  this  means  that 


(x-p^Cx-Pj)  <  0 


(D-81) 


or  x  must  lie  within  the  range  of  the  roots  p^  and  p2«  This  means  that  since 
0  £  *  <_  u,  the  second  root  p2  cannot  be  positive,  because  then  there  will  be  an 
interval  (0^)  for  which  a  value  of  x  will  yield  (x-p^)(x-P2)  non-negative. 


(D-82) 


means 


2 

u 


(D-83) 


At  x  ■  0,  where  the  supports  are,  m^  becomes 


m;c°) 


(D-84) 


For  the  function  not  to  become  lees  than  -1,  m^(o)  must  be  non-negative, 


D-17 
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i«6« ) 


mx(o)  *  I  -  ^  -  0  (D"85) 

or 


2 

u 


< 


6_ 

2 

c 


Since  the  shell  is  short, 


c  <  6  means  that 


(D-86) 


2 

u 


>  1 


(D-87  ) 


However,  at  the  beginning  and  time  t  ■*  1  the  initial  value  that  u  assumes  is  uq 
given  by  Equation  (D-71) ,  i.e., 


2 

u 

o 


_ 6 _ 

c2(p-l) 


which  cannot  exceed  1,  i.e., 


(D-88) 


%  *  ~~2~~ -  -1  or  Pi1  +  1 

c  (p-1)  *  c 

As  time  goes  on,  we  know  that  u  increases,  as  per 

u2  ■  41  +  (M9) 

D-18 
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until  at  t  < 


c2 

“2  P* 

c+6 


u  becomes  1. 


2  6 

In  summary,  therefore,  c  <  6  and  p  >  1  +  -^  guarantees  that  m'(o)  >_  o. 

c 

At  the  other  end,  x  ■  u,xm'(u)  ■  o,  and  it  is  the  second  derivative  that  dictates 
that  the  moment  resultant  is  less  than  or  equal  to  1,  since 


m»(u)«  -  ^  -  c2  -  -  +  c2^  <  0  (D-90) 


Returning  to  the  velocity  equation  we  have,  noting  that 


2  _  6t 
u  -  —x - 

c  (p-t) 


(D-91) 


T  “ 


(D-92) 


(P-t)2 


c  u 


(D-93) 


and 


6  + 


2  2 
c  u 


J2 

P-T 


(D-94) 


D-19 
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D-20 
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D-21 
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(D-102) 


,  (  4 

6  )  uc 

Sb  c4(p-l)2 

2  ,/ti2  2)2 
c  |  2  \6+c  u  1 

c/Fl  72p2 

_  Sb  c4(p-l) 

77)  '  c/Fi  48*> 


(D-103) 


■  _  & 

+  xp 


„  /  ,>3/2 

c  (p-1; 

2/b  p2 


/6  c  /y1! 

8  p 


+ 


(D-104) 


n-29 
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This  can  also  be  rewritten  In  terms  of  time  t  as 


w(x,t) 


f(x)  +  x 


— /7(p -x)3^2  +  ^£-p/x(p-x)  + 

2/6  4/5" 


3c  2 

- p 

k/6 


f(x)  +  X 


2/S”( 


'/t(P-t)3/2  -  (p-l)3/2U 


2£-p  j/TTPH  -  /^|  +  ^P2  j  tarT1^  - 


(D-105) 


Using  the  expression  for  w  In  terms  of  (x,t),  we  observe  that  at  one  end  at 
time  x  ■  1  for  all  x  £  uq,  the  expression  In  brackets  vanishes,  while  f(x) 
assumes  the  form 


f(x)  -  £<p-i>(;r) 
By  Equation  (D-34) 


for 


for 


x  <  u 
—  o 


X  <  u 
—  0 


(D-106) 


(D-107) 


D-23 
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To  obtain  the  solution  for  uq  <_  x  <_  u  we  proceed  as  follows.  At  this  point  we 
need  to  use  the  value  of  w  of  x  ■  u,  which  has  not  been  obtained.  Using  this 
quantity  [(see  second  interval  for  u  <  x  <  1)  (Equation  (D-130)]  we  must  have 


w(x,t)]  -  '2'[2Pt  “  t2 

U“X  L 


-  2px  +  T 


2* 


(D-108) 


Rewriting  all  t's  in  terms  of  u  [Equation  (D-92)]. 


w(x,t) 


U“X 


-  -  £ 
2 


1  -  2p 


UV)  P  UV 


(D-109) 


and  using  Equations  (D-103),  (D-108),  and  (D-109)  we  obtain  Equation  (D-110). 


f(u)  +  up 


2 


/6  c^(p-l)2 

c/fT  72p2 


6  j  uc4 _ /6  c4(p-l)  1 

c2  (8(6+c2u2)  ’  c/FT  48p  ) 


tan 


-1 


2  2 


c  u 


+  p- 


A  A 
c  u 


<64c2u2)2 


(D-110) 


D-2A 
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and 


w(x,x)  -  -  -| 


&  2 
-a-  cp  x 


2  2  4  4 

C  X  .  C  X 

p  '  2  2\2 


1  ‘  2P  (6+c2x2)  +  P  (6+c2x2): 


tan 


■‘(s)  ■  •■"(») 


xp 


3c2u 


,2  ,2 

3c  x  3c  u 


3c2x 


(6+c2u2)2  (6+c2x2)2  4 (6+c2u2)  (6+c2x2) 


for  u  <  x  <  u.  (D-113) 

o  —  — 

This  agrees  with  the  two  end  values.  Furthermore,  w  >  0  for  all  p  >  t  and  w”  *  0 
as  required. 


2.  Second  interval  (u  <  x  <  1)  and  1  <  t  <  t’,  where  x*  is  defined  later. 


Solve 


— ~  m"  +  p  +  n.  -  w  ■  0  (D-114) 

2c2  *  * 


with 


p  =*  0 


(D-115) 


m  ~ 
x 


-1 


(D-116) 


D-26 
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(D-117) 


w  >  0 


(D-118) 


w"  <  0 


(D-119) 


The  equilibrium  equation  reduces  to 


w  “  -1 


(D-120) 


By  integrating  with  respect  to  time 


w  -  -t  +  C, 


(D-121) 


T  -  1  ,  W(l)  »  p  -  1 


CD-122) 


hence, 


p  -  1  -  -1  +  C, 


C1  "  p 


i  •  e  • , 


w(t)  ■  p  -  t 


for  all  p  -  r  >  0. 


(D-123) 


(D-124) 


(D-125) 


(D-126) 
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Integrating  again 


w(t) 


1  2  .  „ 
Pt  -  2  T  +  c2 


(D-127) 


However,  at  t  *  1  by  Equation  (D-55) 


w(l)  -  “(p-1)  *  p  -  ~  +  C2 


(D-128) 


(D-129) 


Hence , 


w(t) 


(D-130) 


STAGE  3  (t’  <  t  <  x  ) 
—  —  o 


We  observe  that  when  u  reaches  the  midpoint  (u  ■  1) ,  the  velocity  w  is 
still  not  zero,  i.e.,  for  time  t'  when 


u 


2 


6t’ 

C2(p-T* ) 


1 


i.e. , 


(D-131) 


(D-132) 


D-28 
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w(t»)  -  (p-r’)x  +  o 

the  hinge  has  reached  the  center  and  motion  continues. 

Again,  we  must  consider  the  possible  portion  on  the  yield  surface, 
concludes  that  portion  AB,  for  which 

a*-'1 

-1  _<  jC  1 

w"  -  0 


(D-133) 


Hodge 


(D-134) 


(D-135) 


(D-136) 


w  >_  0  (D-137) 

is  the  correct  location.  In  addition  to  the  above  we  must  have  continuity  of 
velocities,  accelerations,  displacements,  and  stress  resultants,  at  time  t  *  x'. 

The  equilibrium  equation  to  be  solved  then  is 

— »  nT+p+n.  -w*0  (D-138) 

2c  *  ♦ 

with 

n.  -  -1  ,  p  -  0  ,  i.e., 


D-29 
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—k  a"  ■  1  +  w  (D-139) 

2c  x 


However,  by  Equation  (D-66)  or  Table  D-3  we  know  that  at  t  ■  t* 


w(t') 


t-t  ' 
u-1 


and 


(D-140) 


(D-141) 


Since  the  velocity  must  be  non-negative  for  the  strain  vector  to  satisfy  the 
flow  rule  on  the  yield  surface  and  since  its  second  derivative  with  respect  to 
the  space  variable  x  is  zero,  the  velocity  can  only  be  linear  in  z.  Assume 


w(t)  -  f A1  +  A2t)X 


(D-142) 


then  w(t)  -  A2x  (D-143) 


D-30 


From  the  boundary  condition 


A 


2 


(D-144) 


Hence, 


w(t) 


(D-145) 


(D-146) 


From  the  boundary  condition  on  velocity,  however, 


3 

2 


x 


or 


A 


1 


P 


Hence, 


(D-147 ) 


(D-148) 


(D-149) 


This  equation  is  valid  for  all  times  for  which  w(t)  j>  o,  i.e., 


NSWC  TR  86-328 


P  " 


0 


(D-150) 


i.e.,  for  times  t  In  excess  of  x' 


but  less  than  or  equal  to  tq 


2 

t  <  P  (D-151) 

°  c  +2 

tq  represents  the  time  for  which  velocity  w(t)  vanishes  and  motion  stops. 

Before  we  Integrate  the  velocity  to  obtain  the  displacement  distribution  we 
use.  Integrate  the  equation  of  equilibrium 


_i_  -  i  -  3  c  +21 

,*'■  1  He2'1 


(D-152) 


1  .  -  3/c  +21  2  _ 


(D-153) 


_1_  ,  xf  1  (cjrt)  3  c  - 

2  mx  2  4  £  x  +  Cl  x  +  C2 

2c  c 


(D-154) 


At 


x  ■  1,  m'  ■  o,  i.e. , 


,  _  6^c_ 

'1  /  2 
4c 


(D-155) 


D-32 
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(D-156) 


(D-157) 


m^x-l)  -  1 

giving  C2  ■  -1.  Therefore, 


(D-158) 


(D-159) 


(D-160) 


We  also  note  that  m^Co)  •  -1,  as  It  should. 

Furthermore,  as  Equation  (B-88)  Indicates,  the  yield  criterion  Is  not 
violated.  [See  Appendix  B,  Equation  (B-81),  short  shells,  low  loading.] 


Integrating  the  velocity  w(t)  and  applying  the  boundary  conditions  that, 
in  the  two  regimes  (and  for  t'  £  t  £  tq) 


(a)  0  <  x  <  u  ,  and 

—  —  o 


D-3j 
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(b)  u  <  x  <  1 
o  — 


displacements  at  t  *  must  match,  we  have 


w(t)  =  -|x 


+  E^x) 


3  2 
3*P 


(c2-H0)c2 

(c2+6)2 


+  E^(x) 


(D-161) 


In  the  interval  0  £  x  <.  u0 


at  t  *■  t1 


2 

c  +  6 


(D-162) 


we  have 


»(l,T')  -  i(p-l)  J-  + 


/r7  (p-T1 


(p-l)3/2  \  +  —  p  j  /t'(p-t’)  -  /p  -  1  1  + 

’  4/5"  1 


(D-163) 


Observe  that  at  that  time 


_jJ_  ,  £  (D-164) 

p-x  ’  6 


D-34 
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I 

I 


Hence  for  o  <  x  <  u 
—  —  o 


Ei(x) 


L  ,,x  3  2 

"  2(p_1)r  "  4 xp 

o 


+ 


x 


2 

P 


2  2 

c  /  ,  ,3/2  ,  3c  p  3c  — n  . 

-  (p-1)  +  ■  fo~\ - P  /P“1  + 

2/6  4  (c  +6)  4/5 


For  u  <  x  <  u  and  t 


I 


w(x,t’)  “  -  ^ 


2  2  4  4 

,  c  x  .  _  c  x 

1  -  2p  -t - o~ oT  +  P 


"  (6+c2i2): 


4  cp2  !  j  -  taa-^yj ij 


xp 


-  2  ,  ,.3/2 

3c  *  -  (.P-.1) _  j  + 

(6+c2x2)2  2/5  p2 


■3  2 
3c  x 


4  (6+c2x2) 


A  /FI' 


8  c  p 


xp 


i  2  ,  ,,3/2 

3c  u  c  (p-1) 

2 


y(6+c2u2)2  2/6  p2  l  4 (6+c2u2) 


+  ,  3c  u  -  _  A  c  ^FI 

1  ~  8  p 


_/6  2 
8 


r1/^  _  rar*_ M  ZUH 
(/6 


cp  x  j tan  j—  }  -  tan  f  Ij 


(D-169) 


(D-170) 


D-37 
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or 


(D-173) 


(D-174) 
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Hence, 


We  can 
obtain 


in  the  Interval  u  <  x  <  1 
o  —  — 


E-^x)  - 


3  _ 2  (c2+lo)c2 

'T,p 


£ 

2 


2  2  4  4 

-  n._  C  X  C  X 

1  ”  2p  ,  0  o\  +  P 


U+c2x2)  '  (6+c2x2)2 


A  2 

“jcpx 


•  ““'IS 


Si 


*p 


j-j£ 


3c2x 


.  3c2  3c2x  ) 

(6+c2)2'  (6+c2x2)2  I(^?I'4(6+c2x2)i 


(D-175) 


calculate  the  displacement  at  rest  (x  *  xq)  by  replacing  x  by  xq.  We 


p2x  +  E1 


(x) 


X"X 


(D-176) 
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TABLE  D-1.  SUMMARY,  SHORT  SHELLS.  HIGH  LOADING,  1 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

CONDITIONS 

“o  “  C2(P-1)  0  <r  <  1  0  <x  <u0 

MOMENT 

RESULTANT 

m  <*.r)  =  2[*_-1 

[uo  J 

3t’OR=2^]3-6[e-]2*6[^-’ 

i 

MEMBRANE 

RESULTANT 

DISPLACEMENT 

uo 

VELOCITY 

w  =  (p  ~  1)  T  „ 

U0 

ACCELERATION 

W  =  —  1)  X 

uo 

TIMEt 

o 

-  -  -  - - 

N/A 

DISPLACEMENT 
AT  REST 

HAS  NOT  COME  TO  REST  YET 

D-41 
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TABLE  D-2.  SUMMARY,  SHORT  SHELLS,  HIGH  LOADING,  2 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

SHORT  (c2  v<  6 ij 

■■oai 

CONDITIONS 

o  6  0  <  T  <  1 

uo  =  c2(p-l)  uc<x<1 

MOMENT 

RESULTANT 

mx  =  1 

MEMBRANE 

RESULTANT 

•V  =  -1 

DISPLACEMENT 

w  =  l(p-l)r2 

VELOCITY 

1 

£ 

II 

5 

ACCELERATION 

w  =  p  —  1 

TIME  T 

o 

N/A 

DISPLACEMENT 

AT  REST 

. 

HAS  NOT  COME  TO  REST  YET 

TABLE  D-3.  USEFUL  FUNCTIONS  AND  THEIR  DERIVATIVES 


96  '/6 
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TABLE  D-4.  SUMMARY.  SHORT  SHELLS,  HIGH  LOADING,  3 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

SHORT  jc2  <6^ 

CONDITIONS 

j  6  r  0<x<u  PP 

“  'c2(p-j  1<T<f. 

MOMENT 

RESULTANT 

MEMBRANE 

RESULTANT 

n„  =  -1 

9  6 
‘c2(p-i) 

D-44 
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TABLE  D-5.  SUMMARY;  SHORT  SHELLS.  HIGH  LOADING,  4 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

SHORT  C2  <6 

HIGH  LOAD  p>  1 

C* 

CONDITIONS 

u  <  x  <  1 

2  St  ^  T*  _  _r_E P  =  -E_ 

ts- 

MOMENT 

RESULTANT 

mx  =  -1 

MEMBRANE 

RESULTANT 

DISPLACEMENT 

w|t)  =  JL  j^p  t  -  r2-pj 

VELOCITY 

w  (t  )  =  p  —  r 

ACCELERATION 

w  (  t)  =  -  1 

TIME  T 

o 

AT  TIME  t'  =-E_,  TRAVELING  HINGE  u  MOVES  TO  MIDLENGTH  x  =  u  -  1  AND 

Pi 

THIS  REGIME  SHRINKS  TO  ZERO 

DISPLACEMENT 
AT  REST 

HAS  NOT  COME  TO  REST  YET 
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TABLE  D-6.  SUMMARY,  SHORT  SHELLS,  HIGH  LOADING,  5 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

SHORT  C2  <  6 

HIGH  LOAD  p>  1  +  -L 

C2 

CONDITIONS 

r2 

p  >r  r'  <  t  <  tq  t'  =  °  p 

C2  +  6 

MOMENT 

RESULTANT 

"■  ■  - E'’?]"*  'f-f]"’ 

MEMBRANE 

RESULTANT 

■V"- 1 

DISPLACEMENT 

w(t)  =-^-xfpT  1  (?2-  2Jr2l  +  E-.  (x)  ‘ 

2  L  2  c2  J  1 

WHERE  FOR  0  <  x  <  uQ  E,  IS  GIVEN  IN  TABLE  D-7 
uQ<  x  <1  Ef  IS  GIVEN  IN  TABLE  D-8 

VELOCITY 

"■(2Hp~i£^T] 

ACCELERATION 

TIME  T 

0 

DISPLACEMENT 

AT  REST 

WHERE  E1  WILL  BE  CALCULATED  EITHER 

T  _  _  BY  TABLE  D-7  OR  TABLE  D-8 

'  0 
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TABLE  D-7.  E^x)  FOR  0  <  x  <  uQ 


E  i  ( x  ) 


aJ.(p-1)  JL-3  xpZ  jc^jc2  + 

2  uo  4  (C2  +  6)2 


X 


3C2  p2 
4(C2  +  6) 


-J£_  pv^T  +  ^LcP2 

4n/6  8 
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CASE  D 


APPENDIX  E 


-  LONG  SHELLS,  HIGH  LO i 


2 

c  >  6 


P 


>  I  + 


E-l 
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LONG  SHELLS,  HIGH  LOADING,  (c2  >  6  ,  p  >  1  +  ^ 

This  case  is  summarized  in  Tables  E-l  through  E-7 . 

Let  us  examine  first  the  time  interval  0  <_  t  <_  1.  Table  E-l  (it  applies  for 

0  <  x  <  u  )  and  Table  E-2  (it  applies  for  u  <  x  <  1)  summarize  the  results. 

—  —  o  o  —  — 

The  pertinent  intervals  on  the  yield  surface  lie  along  AB  and  at  point  B 
(corner),  respectively.  The  analysis  for  "short  shells  under  high  loading," 
therefore,  applies  for  times  0  <  x  <  1.  Equations  (D-36)  and  (D-37)  will  be 
satisfied  and  hence  there  will  be  no  violation  of  the  yield  locus.  Furthermore, 
the  velocity  distribution  from  Equations  (D-29)  and  (D-20)  is  such  that  the  flow 
rule  is  satisfied.  Therefore, 


w(x,t) 


(p-1) 


0  <  x  <  u 
—  —  0 


(E-l) 


w(x,t)  *  (p-l)t 


U  <  X  <  1 
o  —  — 


(E-2) 


At  time  t  *  1  the  load  is  removed.  We  observe  that  the  analysis  of  the  previous 

2 

section  for  the  next  time  interval  (t  1)  is  only  valid  for  c  <  6.  For  long 
2 

shells,  however,  c  >6. 
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Hence  a  new  assumption  is  required  as  to  the  ranges  on  the  yield  surface. 
This  must  be  done  in  association  with  the  fact  that  in  the  interval  o  £  t  <  1  the 
velocity  profile  is  given  on  Figure  E-l  (this  figure  comes  from  Figure  D-l  appli¬ 
cable  to  short  shells  and  high  loads).  Figure  E-l  suggests  that  for  1  <_  r  <  tq 
the  logical  compatible  velocity  profile  is  given  by  Figure  E-2  in  such  a  way  that 

1.  at  time  t  ■  1,  y  -  0 


2.  at  time  t  -  l,  u  •  u  and  the  range  AD  has  disappeared. 

Based  on  the  previous  assumptions  we  must  solve  the  equilibrium  equation  in 
three  intervals  for  1  £  r  <  tq,  and  account  for  Initial  and  boundary  conditions. 

1.  For  times  1  £  t  _< 

a.  For  0  £  x  £  y  (Range  AD,  Figure  A-3,  Table  A-l) 


w  ■  0 


w"  >  0 


+  p  -  w 


0 


p  «  o 


(E-3) 

(E-4) 

(E-5) 

(E-6) 

(E-7) 
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b.  For  y  £  x  <  u  (Range  AB,  Figure  A-3,  Table  A-l) 

%  “  _1 
v  >_  0 

w"  =  0 

— +  n.  +  p-  w=»0 
2c  X  * 


p  *  0 


c.  For  u  £  x  £  1  (Point  B,  Figure  A-3,  Table  A-l) 


m  =  1 
x 

v<  >  0 


w”  <  0 


— rin”  +n+p-w*0 
„  z  x  $ 


p  *  0 


(E-8) 

(E-9) 

(E-10) 

(E-ll) 

(E-12) 

(E-13) 

(E-14) 

(E-15) 

(E-16) 

(E-17) 

(E-18) 
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We  know  that  at  time  r  ■  1,  when  the  pressure  load  ceases  to  act  (see 
Tables  D-5  and  D-6) ,  the  Initial  conditions  are 


1.  For 


<  x  <  u  ■  /- 

-  -  °  v« 


c2(p-l) 


w(x,l)  - 


w(x,l)  -  | 

o 


2.  For  u  <  x  <  1 
o  —  — 


w(x,l)  -  p-1 

w(x,l)  ■  ~  (p-1) 

We  now  examine  the  time  Interval  t  1. 
three  segments 

1.  0  £  x  <_  y, 

2.  y  <  x  <  u,  and 


3.  u  <  x  <  1 


(E-18a) 


(E-18b) 


(E-19 a) 


(E-19b) 


The  space  Interval  is  subdivided  in 


Both  y  and  u  will  be  defined  later 


E-5 
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FIRST  INTERVAL  0  <  x  <  y 


The  equilibrium  equation  yields 


n,  “  w 


Since  w  =  0 


for  all  times  t  >  1,  however,  this  means  that 


w  -  0 


and,  hence. 


%  "  0 


and  integrating  Equation  (E-21)  once  with  time 


w  *  C^Cx) 


Hence, 


1  (p-1) 

w  =  2  u - x 

o 


In  summary,  for  0  _<  x  _<  y 


m  ■  -1 
x 


E-6 


(E-20) 


(E-21) 


(E-22) 


(E-23) 


(E-24) 


(E-25) 


(E-26) 
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n.  ■  0 
9 


if  -  0 


w  *  0 


w 


i  <£li>x 


2  u 


and  y  has  not  been  determined. 


(E-27) 

(E-28) 

(E-29) 


(E-30)  ] 

i 


SECOND  INTERVAL  y  <  x  <  u 

j 

I 

The  equilibrium  equation  yields 

(E-31)  j 


Since  by  Equation  (E-10)  w"  -  0,  the  velocity  profile  must  be  linear  in  x. 


Assume 


»  x  +  w 

2x  x  t  w 

2c  X 


w  -  A(t)x  +  B(t) 


(E-32) 


Hence , 


w  ■  A(t)x  +  B(t) 


(E-33; 
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Hence, 


in"  -  1  +  A(x)x  +  B(x) 
2c  x 


Integrating  once  with  respect  to  x 


-in'  =  ^-iT— x  +  (B(t)+1)x  +  D(x) 
2c  x  L 


and  once  more 


-i^m  “  ^A(t)x3  +  i(l+B(x))x2  +  D(t)x  +  E(t) 
2c  X  °  Z 


At  x  =  y  the  moment  resultant  and  shearing  forces  must  agree,  i.e., 


-i^m'(y,x)  -  y2  +  (l+B(-r))y  +  D(t)  "  0 
2c  L 


(y,x)  ■  ^A(x)y3  +  i(l+B(x))y2  +  D(x)y  +  E(x)  - - 

2c  X  0  L  2c 


Also  at  x  ■  u 


— ^m'(u,x)  -  4a(x)u2  +  (1+B(x))u  +  D(x )  -  0 
2c 


E-8 


(E-34) 


(E-35) 


(E-36) 


(E-37) 


(E-38) 


(E-39) 
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(u,t)  *  -tA(t)u3  +  4(1+B(t))u^  +  D(t)u  +  E(t)  =  - 

2c2  X  6  2  2c2 


(E-40) 


Solving  Equation  (E-37)  for  2c  D,  we  get 


2c2D  -  -c2A  y2  -  2c2(l+B)y 


(E-41) 


2  2 

Solving  Equation  (E-38)  for  2c  E  and  replacing  2c  D  from  Equation  (E-41),  we 
obtain 


2c^E  =  -1  +  |c2Ay3  +  c2(l+B)y 


(E-42) 


Replacing  2c  D  in  Equation  (E-39)  from  Equation  (E-41)  we  obtain 


:2X(u2-y2)  +  2c2(l+B)(u-y)  =  0 


(E-43) 


2  2 

Replacing  2c  D  and  2c  E  in  Equation  (E-40)  from  Equations  (E-41)  and  (E-42), 
we  obtain  a  second  relation 


ic2A [u3-3uy2+2y3]  +  c2(l+B)(u-y)2  -  2 


(E-44) 


^c2A(u-y) (u2+yu-2y2)  +  c2(l+B)(u-y)2  =  2 


(E-45) 
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Substituting  c^(l+B)(u-y)  from  Equation  (E-43)  in  Equation  (E-45)  we  finally 
obtain  the  value  of  A,  i.e., 


A 


12  1 

2  f  v3 
c  (u-y) 


(E-46) 


Hence  (for  u  ^  y) 


2(1+B)  -  -A(u+y) 


which,  in  view  of  Equation  (E-46),  becomes 


2 (1+B)  -  -(u+y)x/^\— — 
\C  / (u-y) 


or 


1+B 


6_  (u+y) 

2  ,  v3 
c  (u-y) 


Therefore, 


12  1 

2  ,  ,3 

c  (u-y) 


1+B 


6_  (u+y) 

2  ,  x3 
c  (u-y) 


» 

I 


(E-47) 


(E-48) 


(E-49) 
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B  -  K,  -1 

c  (u-y) 


2c2D  -  ^ 
(u-y)' 


(y+u) (y^-4yu+u2) 
Cu-y)3 


Also  by  Equation  (E-41) 


2c2D - 


(u-y)’ 


By  Equation  (E-42) 


2c2E  -  -1  4-  22-13M1 
(u-y) J 


m  (x,t)  =  — - — ^ f-4x3+6(u+y)x2-12uyx-(y+u) (y2-4yu+u2)] 
1  (u-y)3L 


m'(x,x)  =  — - — J-12x2+12(u+y)x-12uy 
X  (u-y) J  L 


m"(x,x)  ”  - r[-2x+(u+y)] 

x  (u-y)3 


(E-50) 


(E-51) 


(E-52) 


(E-53) 


(E-54) 


(E-55) 


(E-56) 


(E-37 ) 
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We  observe,  that,  in  fact 


OjCy.x)  +  -1 
m^(y,t)  *  0 

n^Cujx)  1 
m^Cu.r)  -«■  0 


■J(y.T) 


12 

(u-y)2 


m^(u,T) 


12 

(u-y)2 


Before  we  can  combine  all  three  ranges,  we  examine  the  third  interval  (u  <  x  <  1). 


THIRD  INTERVAL  (u  <  x  <  1) 


The  equilibrium  equation  yields 

w  -  -1  (E-58) 


Integrating  once  we  get 

w(t)  -  -t  +  Cx  (E-59) 


E-12 
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However,  by  the  initial  condition  in  velocity  [see  Equation  (E— 18) ]  at  t  *  1, 
we  have 


[-t  +  Ci]  -  p  -  1 

T-l 


(E-60 


Hence, 


Cj_  “  P 


and  w(x,x)  “  p  -  t 


Integrating  again 


(E-61 


(E-62 


w(x,x) 


1  2 

px  -  -t  +  D1 


(E-63 


By  Equation  (E-19)  at  t  ■  1 


(E-64 


Leaving  Equation  (E-63)  with  the  known  constant  we  observe  that  the  displace¬ 
ment  must  agree  with  that  one  obtained  from  the  second  interval  for  all  times 
r  >  1.  We  must  also  retrieve  the  value  of  the  displacement  at  time  t  *  1  as 
given  by  Equation  (E-19),  i.e., 


w(  x ,  1)  “  -|(p-l) 


(E-65 


E-13 
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and,  therefore,  takes  on  the  value  already  given  by  Equation  (E-64)  and 

w(x,t)  “  p(t  "  7)  "  (E-66) 

It  can  easily  be  seen  that  for  the  time  range  of  Interest,  w(x,t)  Is 
non-negative. 


Vfe  proceed  to  match  solutions  in  the  three  intervals  for  t  >  1, 

since 

1. 

0  £  x  <  y  w-0 

(E-67) 

2. 

y  <  x  <  u  it  *  A(t)x  +  B(t) 

(E-68) 

3. 

u<x<1  w*p-r 

CE-69) 

For 

x  ■  y  A(t)y  +  B(t)  -  0 

(E-70) 

For 

x  ■  u  A(t)u  +  B(t)  *  p  -  t 

(E-71) 

Therefore, 

B(t)  •  — A(t )y 

(E-72) 

A(T)(u-y)  *  p  -  x 

(E-73) 

(E-74) 
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Differ satiate  Equation  (E-74)  with  time 


A(t) 


'  Tu^yT 


(p-t) 

(u-y) 


2  rr'“-y> 


(E-75) 


By  Equation  (E-48)  however  and  setting  9  ”  u-y 


A 


12  1 
c  (u-y) 


i.e. , 


(E-76) 


we  get 


1  (p"T>  d  „  12  1 

'  9  “  “2“  JF6  'll 

9  C  0 


(E-77) 


or 


(P-t)  d_fl  _  1  12  1_ 

a2  dx”  0  "  2  .3 

0  c  0 


(E-78) 


(E-79) 


(E-80) 


1  [~12  1 
(p-O  c2  (u-y) 


(E-81) 


E-15 


NSWC  TR  86-328 


Taking  the  positive  root  only  since  u  >  y  we  obtain 
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Equation  (E-50)  gives  B 


3  >  L  (u+y>  _  i 

D  0  'X  x 

c  (u-y)J 


Also  Equation  (E-71)  relates  A(t)  to  B(t) 


A(t)u  +  B(t)  =  p  -  t 


(E-95) 


Differentiating  (E-71)  with  respect  to  time  x  we  obtain 


A(t)u  +  A(t)^  +  B  =  -1 


(E-96) 


By  Equation  CE— 74) 


A(x  ) 


(p-t) 


(E-96a) 


By  Equation  (E-48) 


A(t )  =  - 


12 

2 


(u-y)3 


(E-96b) 


Replace  A(t)  from  Equation  (E-48),  A(t)  from  Equation  (E-74),  and  B  from 
Equation  (E-50)  in  the  differential  Equation  (E-74) 


12  u  (p-r)  du  6  (u+y) 

2  .  .3  dt  2  ,  ^3 

c  (u-y)  }  c  (u-y) 


(E-97 ) 
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or 


! 


(p~t)  du  _  6  1 

iT  c2  (u-y)2 


(E-98) 


or 


du  *•  — 


/6 

dx 

c 

*  “ 

1 

(p-t) 

2  .  ^(p-o2 

2 

(p-l)J 

dx 

“  - 

1 

,  ^  /§ 
(p_T) 

2  -  «E^(p-t)2 

2 

C 

(P-1)3 

and  by  Equation  (E-93) 


du 


6  dx 

^  TF^Te 


CE-99) 


By  Equation  (E-93)  6  is  given  In  terms  of  x 


9 


(^<p-x)2 

(p-ir 


1 

2 


de 

dx 


(2p-3) 

(p-1)3 


2 


-  <2£z3)(p-t,2 

(p-l)J 


1 

2 


(E-100) 


E-19 
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de 

dr 


/6  \  (2p-3)  (p-r) 
\c2)(p-l)3  8 


CE-101) 


and  since 


de 


^  Wo 

(p-1)3  c 


dt 

e 


(E-102) 


dt 


ede  (p-l)3  c2 

(p-x)  (2p-3)  6 


(E-103) 


and 


du 


6_  1  1 
c2  (p-x)  e 


ede  (p-1)3 
( P~t)  (2p-3) 


% 


1  (p-1)3 

(p-x)2  *2Fsy 


de 


(E-104) 


2 

and  solving  Equation  (E-94)  for  (p-x)  and  replacing  it  in  Equation  (E-104), 
we  get 


(p-x)2 


(p-1)3 

(2p-3) 


2 


2 

. 


(E-105) 


(P-x)2 


(2p-3) 

(p-1)3 


1 


2 


2 


(E-106) 


du 


(2p-3) 


a  r 


de 


de 


T  n-T- -  on  (E-107) 
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/3 

2c 


de 


+  c 


l 


at  t  ■  1 


u  - 


A, 


(u.  *  ¥) 


(uo  -  ¥) 


Hence , 


Integrating  once  with  respect  to  e  we  get 


u 


(E-107a) 


(E-107b) 


(E-107c) 


(E-107d) 


(E-108) 
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(1_Uo) 


(  ,  2/3\  — - — 

\  o  c  /  /3 

-  e 


I  2/3 

u - 

\  o  c 


1  c 


2c/l-u  \ 

2/3\  -i - 21 

Uo  c  /  /3 

-  e 


(“o  -  ¥) 


Multiply  top  and  bottom  with  e  uo^,/^  to  obtain 


Lu  cosh 
° 


u  sinh 
o 


m :  f  -“(^l 

m  *  f 


For  the  general  case  of  any  e,  using  Equation  (E-110)  we  obtain 


L  co6hpj: 


si)  + 


(U~U.>)\  |  2/3 


"-"o)\ 

A  /_ 


c  /c  (u-u 

u  sinhl— — — 
°  \  /3 


Equation  (E-120)  can  be  written  in  the  alternative  form 


LJMV  2/5 

2/S  .  \  A  )  "V 


(E-119) 


(E-120) 


(E-121) 


(E-122) 
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Solving  Equation  (E-100)  for  we  obtain 


T1 


P 


A 


(p-1); 


(2p-3) 


(E-123) 


We  distinguish  three  cases 


2 

c  2  o  o 

1*  2  -  *  o  (it  will  turn  out  that  p  ■  j  and  c  >  12) 


(E-124) 


2.  2  - 


rsi >  0 


for  p  >  j 


(E-125) 


3.  2  - 


c2  2 

g-0  <  o 


for  p  <  j- 


(E-126) 
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CASE  1 


We  now  examine  case  6 


2/3 


The  differential  Equation  (E-82)  reduces  to 


-  (p-T)efi 


(E-127) 


or 


de 

dr 


Since  p  i1  t  and  e  +  0 


Equation  (E-98)  reduces  to 


,  6  1  dr 

U  c2  (p-r)  0 


when  0  is  a  constant,  i.e., 


(E-128) 


(E-129) 


du 


6_  J_  dt  _  /3  dr 

c2  2/3  c  ^P_T^ 

c 


(E-130) 


E-26 


NSWC  TR  86-328 


and  Integrating 


u  - loge|p~r|  +  logeC 


At  time  t  ■  1,  u  •  uQ,  i.e.. 


A 


o  c 


loge  p-1 |  +  logeC 


} 


A 


lo*elt^n 


+  u 


And 


y  “  u  -  6  =  u  -  0  + 
1  o 


A  (p- 
c  °8e  (p-1) 


=I>I  - 


M  +  log 

c/p-1  c  c  e 


<2l4 

(p-iT 


However,  at  r  ■  1,  y  *  0.  This  leads  to  p  *  y  and  since 


P  >  1  +  ^ 

c 


c  >  12 


(E-131) 


(E-132) 


(E-133) 


(E-134) 


(E-135) 


1 


1 

*1 
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By  Equation  (E-74) 


A(t) 


(p-t) 

u-y 


(p-x) 

e 


c(p-t) 

2/5 


(E-136) 


By  Equation  (E-48) 


X(r) 


12  1  12  1  c 

c2  (u-y)3  c2  8 3  2/5 


(E-137) 


By  Equation  (D-72) 


B(t)  =*  -  A(x)y  ”  -  —  P~-T^-y 

2/3 


(E-138) 


By  Equation  (E-50) 


b  -  L.  (u+y>  _  i  ,  _S_ 

c2  (u-y)3  4/5 


O  o  ^  2/5. 

2u  -  8  +  - log 

o  c  e 


EH 

P-1 


-  1 


(E-139) 


Therefore, 


P 


3 

2’ 


p  >  T, 


c2> 


12 


(E-140) 


u 


o 


/6 


c 


(E-141) 


8 


(E-142) 
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t 


»  ■  “o  +  4l08.|{p=T7 


8  ■  “o  -  9  +  4l0*elT^I7 


(E-143) 


(E-144) 


and  adding  and  subtracting  e  Inside  the  square  brackets  &  can  be  written  as 


a  -  _c  _  1 
B  r?  2 

2/5  1 


(E-145) 


Since 


B(t)  - — (p—r)y 

2/5 


(E-146) 


‘  • '  h +  (p-T)H 


2/5 


y  -  2 


(E-147) 


as  above. 

Also  by  Equation  (E-68)  the  velocity  profile  In  the  second  interval 
y  <_  x  <  u  becomes 


w  ■  A(x)x  +  &(t)  “  -^-(p~r)(x-y) 

2/5 


(E-148) 


and  the  acceleration 


“  X(t)x  +  B(t)  - - —x  +  -£-y  -  i  ■  — “(-x  +  y)  -  -r  (E-149) 

2/5  2/5  2/5  z 
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When  u  reaches  the  midpoint,  let  the  time  taken  be  indicated  by  Then 

we  will  have 


u 


1 


u 

o 


(P-Tl) 

(p-1) 


(E-150) 


1*6* 


i 


(p-Ti) 

(p-i) 


(E-151) 


i«6«  ) 


[p-l]  e 


(E-152) 


If 


P  > 


(E-153) 


p  -  (p-l)  e 


(E-154) 


(E-155) 
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If,  on  the  other  hand,  p  < 


Ti  "  p 


(p-1)  e 


— (1-u  \ 

✓5'  0> 


(E-156) 


p  +  (p-1)  e 


-(1-u  ) 

A(  0> 


p\l  +  e 


(E-157) 


Hence , 


p-%) 

If  p  >  t^,  ■  p  -  (p-1)  e  (E-158) 


If  p  <  r^,  ■  p  +  (p-1)  e  (E-159) 

and  always  ^  1 


It  turns  out  that  the  conditions 


p  >  tx  >  1 


(E-160) 


reduce  to 


p  >  p  -  (p-1)  e 


—  (1-u  ) 

A(  oj 


>  l 


(E-161) 
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or 


0 


which  Is  satisfied 


(E-162) 


and 


1-u 

o 

<  1  (E-163) 

which  cannot  be  satisfied  for  0  <  u  <1 

o  — 

3 

Hence  for  P  =  p  >  t^,  and  >_  1,  the  postulated  profile  violates  the  flow 

3 

rule.  While  for  p  *  j,  p  <  t^,  and  1 


1-u 

o 

>  0  (E-164) 

is  satisfied  for 

0  <  u  <1. 

—  o  — 

Equation  (E-148)  for  the  velocity  profile  and  Equation  (E-9)  (w  >_  o)  imply 
that  the  inequality  is  violated  and,  therefore,  a  different  velocity  profile 
must  be  assumed  for  this  case.  This  particular  case  will  not  be  further  studied 
here,  since  other  considerations,  such  as  inclusion  of  nonlinearities,  have 
priority. 
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CASE  2 


Corresponding  to  a  value  of  0,  call  it  0^  (for  time  t  ■  and 
u  “  If  there  exists  a  value  of  y,  call  it  y^,  such  that 


yi  -  i  -  *i 


(E-165) 


i  •  e  •  , 


O<l-01“y1<l 


(E-166) 


1  >  01  >  0 


(E-167) 


By  Equations  (E-125)  and  (E-167),  therefore, 


0X  <  min  1,  ~ 


(E-168) 


with  p  >  j 
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CASE  3 


Corresponding  to  Equations  (E-126)  and  (E-167)  we  must  have 


2/3 


<  e1  <  l 


CE-169) 


for 


P  <2 


Consequently,  we  combine  cases  2  and  3,  and  treat  them  as  a  single  case,  i.e., 
for 


P  ?  2  an£*  P  >  t. 


We  proceed  now  to  calculate  the  displacement  distribution  for  times  t  such 
that  1  £  t  _<  and  are  given  by  Equation  (E-123).  For  this  purpose,  we 

observe  that  at  time  t  *  1  the  displacement  has  an  initial  value,  depending  on 
the  position  of  the  point  x  at  which  it  is  considered.  This  initial  value  is 
given  by  the  following  table: 
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Therefore 


1.  For  0  <_  x  <  y  and  t  1,  and  by  Equation  (E-21) 


w(x,t)  “  i  -jjp'-x 
o 


2.  For  y  £  *  £  u0  and  t  >_  1 


T 

w(x,t)  -  i  -p--x  +  f  wdT 

0  T-l 


However,  by  Equations  (E-32),  (E-74),  and  (E-72) 


w(x,t)  -  (p-t) 


and,  therefore,  Equation  (E-171)  becomes 


s  .  l  Cp-1)_  .  c2  (p-1)3 

»U,t)  j  +  r  ^7 


<*-uo)  (9-uo) 


lf2'Uo2)  '  -2? (e-“o)  l0*« 


(%-¥) 


(-. ♦  ¥) 


A 

2c 


(9  *  ¥Kl  (9  +  ¥)l  -  (•=  +  ¥KI  (“o  +  ¥)l 


(»-¥KI(‘-¥)M“«-¥KI  (%-¥)!' 


(E-170) 


(E-171) 


(E-172) 


(E-173) 
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3.  For  uQ  <  x  <  u  and  t  1,  the  expression  for  w(x,t)  is  like  the  one 

given  by  Equation  (E-173)  except  that  instead  of  — ^-x ,  the  first  term  is 

1  U° 


4.  For  u  <  x  <  1  and  t  >  1  we  must  proceed  as  follows: 

At  time  t  -  r*,  when  x  -  u  •  u*  and  9  ■  0(u)  “  9(u*)  *  9*  the  displacements  in 
the  third  and  fourth  intervals  must  agree.  This  means  that 


w(u*,t*) 


* 

wdx  + 


(E-174) 


However,  by  Equation  (E-62) 


w(u* ,t* ) 


T*T* 


(p-x)dT  + 


(E-175) 


and  also 
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w(u*,t*) 


¥<-»  +  r  T2FJT  (»*-"„)  (•*-».)  +  Ik2-"2)  - 


(  _  2/3 

/3  /„*  * ,  \o  c 

(0*-u  |  log  - 

2c  l  o)  8e  2/3 

\o  +  c 


-  £  ik  +  ¥KI  (•*♦¥)!- 


:^¥Klk*¥)l>-¥KI  (•*-¥)!♦ 


2/3 

uo  c 


)l08el  ( 


_  2/1 

uo  c 


— 

)|j  =  p(t*-l)  -  j  (x*2-l)  +  Cx 


(E-176) 


Equation  (E-176)  determines  in  terms  of  0*,  u*,  and  x*. 


Therefore, 


w(x,x)  *  j(p-l)  +  p(x-x*)  -  i(x2-x*2)  + 


(“*-%)  (9*-%)  ¥(°*2-uo) - 


k  -  ¥) 


V^(e*-u  )  log  — - =-*■ 

ek  +  ¥) 


-  £  Ik ♦  ¥Ke|  (» 


k +  ¥)i”*el  k +  ¥)l  -  k  -  ¥)io*J  (•*  -  ¥)l + 


“.-¥KI(“o-¥)ll 


(E-177 ) 


and  will  be  given  by  Equation  (E-113) ,  (E-120),  or  (E-122),  while  t-  will  be 
obtained  by  Equation  (E-114)  or  (E-123).  Equation  (E-lll)  defines  y^.  At  that 
instant  in  time,  the  velocity  profile  is  given  by  Figure  E-3.  A  different 
assumption  must  be  made  for  the  motion  to  continue.  This  must  be  so,  since 
y(-r),  where  the  other  hinge  is  located,  varies  with  time  and  has  not  reached  t< 
midpoint  yet.  We  must  also  satisfy  all  geometrical  inequalities,  such  as  y  >  0 
and  0  <  9  <  1  for  the  solution  to  be  valid. 


E-38 


NSWC  TR  86-328 


2.  T-  <  T  <  T 
1  —  —  O 

The  analysis  of  long  shells  under  low  loading  (Appendix  C)  applies  in  this 
case.  We  observe  that  the  initial  conditions  are  different  then  the  ones  in 
Appendix  C.  Both  velocity  and  displacement  profiles  must  agree  for  times  t  ■ 
which  represent  our  starting  time  for  this  interval. 

We  consider  two  intervals. 

a.  For  0  x  <_  y  (along  AD  on  the  yield  surface) 

we  must  have 

w  -  0  (E-181) 

w"  -  0  (E-182) 


m  -  -1  (E-183) 

x 

p  -  0  (E-184) 

and  by  Equation  (E-3) 

*(x,t  )  -  0  (E-185) 
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and  by  Equation  (E-25) 


»<x,n)  -  i 


(E-186) 


Also  the  equilibrium  equation  is 


— ?m"  +  n  +  p-  w**0 
2c2  x  * 


(E-187) 


Therefore, 


%  "  ° 


(E-188) 


w  =  w  =  o 


(E-189 ) 


and,  hence.  displacement  is  independent  of  time  t,  but  dependent  on  location 
x,  i.e.. 


w(x,t)  *  C(x) 


(E-190) 


This  constant  is  the  value  of  the  displacement  from  the  previous  time  range,  i.e. 


CC)  -  i  ifiix 

O 


(E-191) 
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Thus, 


1  <Ei) 


2  u. 


b.  For  y  <  x  <  1  (along  AB  on  yield  surface) 


we  must  have 


w  >  0 


w"  -  0 


n.  «  -1 


-1  <  m  <1 
—  x  — 


p  -  0 


with  initial  condition 


x-y 

(Tl), 

'l-y 

Tl)1 

We  must  satisfy  the  equilibrium  equation 


— rm"  +  n  +  p-  w*0 
2c  x  ♦ 


(E-192) 


(E-193) 

(E-194) 

(E-195) 

(E-196) 

(E-197) 


(E-198) 


E-41 


(E-199) 
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and  assuming  a  distribution  of  velocity  of  the  form 

w(x,r)  *  Ax  +  B 

w(x,t)  ■  Ax  +  B 

and,  hence,  replacing  w  in  Equation  (E-199)  we  get 

-V"  -  1  *  Ax  +  B 
2c2  X 

As  before,  integrating  we  obtain 

-  2c2x  +  2c2[|ax2  +  fix]  +  C]L 
mx  -  c2x2  +  c2|jXx3  +  Bx2]  +  Cx  x  +  C2 

We  must  satisfy  the  following  boundary  conditions 
mx(y,r)  “  -1 
n;(y,T)  “  0 

mx(l,T)  M  1 

m’(l,T)  -  0 
x 


(E-200) 

(E-201) 

(E-202) 


(E-203) 

(E-204) 

(E-205) 

(E-206) 

(E-207 ) 

(E-208) 
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Therefore, 


c2[^&  y3  +  (1+B)y2]  +  CLy  +  C2  -  -1 

(E-209) 

c2[&  y2  +  2(1+B)y]  +  Cx  -  0 

(E-210) 

c2[X  +  2(1+B)]  +  CL  -  0 

(E-211) 

c2[i*  +  ( 1+B )]  ♦  C,  +  C2  -  1 

(E-212) 

After  similar  operations,  we  obtain 

*  ’ '  2<1+B>(^  ‘ '  2<1+4>n^ 

(E-213) 

and  for  y  i1  1 

i  +  b  -  -  -|a+y)A 

(E-214) 

and 

A  -  -  12  1 

c2  (1-y)3 

(E-215) 

1  +  t  -  6,  ^5 
c2  (1-y)3 

(E-216) 

and 

12 

ci  "  -  3y 

1  (1-y)3 

(E-217 ) 
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r  -  i  2(l-3y) 
L2  x  3 

(l-y)J 


(E-218) 


The  bending  moment  distribution.  Equation  (E-204),  becomes 


m  *  1  +  — - — ij[-4x3  +  6(l+y)x2  -  12y  x  -  2(l-3y) 
x  (l-y)3L  J 


(E-219) 


When  x  *  y  (for  t  >  t^)  the  velocity  vanishes,  i.e., 


w(y ,t )  =  0 


(E-220) 


or 


A(t )y  +  B(t )  ■  0 


(E-221) 


However,  y  is  a  function  of  time  t.  Similar  differentiations,  as  with 
Equation  (C-36)  gives 


Ay  +  Ay  +  B  =*  0 


or 


A 


1  _  6_  1 
^  c2  y(l-y)2 


and  the  differential  equation 


(E-222) 


(E-223) 


(E-224) 
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or  y  ■  0  (E-225) 

since  y  +  0 


Therefore,  y(x)  is  a  linear  function  of  time  x  and  such  that  at  time  x  “  x^ 
it  equals  y1  (yCx.^) 


y(x)  -  e1(t_t1)  +  E2 

(E-226) 

E2  ‘  E(,l)  ‘  *1 

(E-227 ) 

yCO  ■  Ex 

(E-228) 

Comparing  the  x  coefficient  of  Equation  (E-198),  which  in  fact  is  A,  with  Equation 
(E-223),  for  times  t  ■  x^,  we  have 


6_ 

2 

c 


or 


But 


(E-229) 


(E-230) 


Ex  -  y  -  yL  (E-231)  j 

i 
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We  also  observe  that  the  velocity  of  the  hinge  circle  is  given  by 


y(f)  =  y1  - 


1 


6  1  I  (>yl) 

c2  (l-yi)2J  P"l) 


(E-232) 


Also  at  time  t 


the  velocities  for  all 


y  <  x  <  1  must  agree,  i.e., 


*(x>ti)  ■  mm*  +  Mti) 


(E-233) 


i.e. , 


A 


(E-234) 


and,  therefore. 


(E-235) 


(E-236) 


and  the  velocity  distribution  becomes 


w  *  A(t)x  +  B( t )  *  A(T)(x-y)  * 


(p-,i) 


(i-y, ) 


6_  (x-y) 
c  (1-y) 


(E-237) 


NSWC  TR  86-328 


Except  the  case  when  either  x  ■  y  or  p  ■  the  velocity  becomes  zero  when 


1 


y 


Thus 


w(x,t)  ■  - =4 


(1-y) 

The  acceleration  is  given  by 


NO  -6] 


„  -  NO 

NO  1  -  *1 TT^ 
L  c  NO 


6  d 

~“c 

l 

c 


■y  "  2  dT 


(x-y) 


(1-y)' 


\  +  6_ 

2 (x-y)  1 

1  2 
c 

(1-y)3  (1-y)2 

•m 

(E-238) 


(E-239) 


For  consistency  w  ^  0.  Since  1  >  y^,  1  >  y,  and  x  >  y,  we  must  also 


have 


P  > 


T 


1 


1  - 


6_ 

2 

c 


(1-y)2 


> 


0 


and 


i  •  e  •  f 


y  <  i 


A 


yi<  1 


A 

c 


(E-240) 


(E-241) 


(E-242) 


(E-243) 

(E-244) 
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At 


y  -  1 


(E-245) 


the  velocity  becomes  zero. 


The  time  taken,  indicated  by  t  ,  is  given  by 


T 

O 


t!  +  - 


(p-Ti)(i-yi) 


.  x  /6 
1_yi  +  ~ 


(E-246) 


Table  E-6  summarizes  the  results  for  the  interval  t.  £  t  £  tq  and  points  lying 
in  0  £  x  £  y. 

The  objective  here  is  to  obtain  the  displacement  distribution  in  the  interval 
y  £  x  £  1  for  times  t  £  x£.  To  do  this  we  observe  that  at  x  *  y  the  displace¬ 
ments  must  be  equal.  Therefore 


w(y,x) 


2  u  y 


(E-247) 


Using  Equation  (E-237)  written  in  a  different  format 


(i-Yl) (p-tj) 


,yi)  ip'Ti)  rr 

l ' 


6 

(l-y)2. 


(x-y)dr 


(E-248) 
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we  get 


Ai  ‘  *-*i 


r,._ 

6  1 

(i-yi) 

L  <2 

C'Ti) 

-  Ti  - 6  1  1  - 

[  =2  h  )2J  i 

I^l) 

-  1  -  H  *  I"1  -  h  ■ 

1 

L  c  1 

(^l) 

(1-yi) 


_  (1-y1)(p-ri) 

D  »  - - — «r 

1 


where  ”  1  -  x 


(E-249) 


(E-250) 


(E-251) 


(E-252) 


(E-253) 


Finally,  Equation  (E-247)  assumes  the  form 


w(y,t)  -  \  ^yK)  +  \  [c2|ai  (t'ti)  -  K(t2‘ti) 


I^loge 


(crBiT)  .  d-x) 
(Wi)  Bi 


(crBiT)  (Wi)  IJ 


(E-254) 


and  the  displacement  at  which  the  velocity  vanishes  is  given  by 


/6 


w/y  ,t  ) .  where  y  “1  - 

Vo*  o/  ’  Jo  c 


and  t 0  is  given  by  Equation  (E-246). 


Table  E-7  summarizes  these  results. 


E-49 


NSWC  TR  86-328 


w 


RANGE  AB  POINT  B  ON  YIELD  SURFACE 

ON  YIELD  SURFACE 


FIGURE  E-1.  VELOCITY  PROFILE  FOR  LONG  SHELLS  AND 
HIGH  PRESSURES  FOR  IN  0  * '  r  <  1 
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POINT  B 
ON  YIELD 
SURFACE 


FIGURE  E-2.  ASSUMED  VELOCITY  PROFILE  FOR  LONG  SHELLS 
AND  HIGH  LOADS  FOR  1  <  r  <  r„ 
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N 


AB  RANGE  ON 


RING 


MIDDLE 


FIGURE  E-3.  VELOCITY  PROFILE  FOR  LONG  SHELLS  AND 


HIGH  PRESSURES 


FOR  1  <  r  <  T, 
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TABLE  E-1.  SUMMARY,  LONG  SHELLS,  HIGH  LOADING.  1 


TYPE 

SHELL  TYPE 

AB  PRESSURE  LOADING  TYPE 

LONG  C*  >  6 

HIGH  LOAD  p  >  1  +  — 

C2 

CONDITIONS 

0  «  T  1 

“°  '  C2(P-1)  0,x<u0 

MOMENT 

RESULTANT 

mx(x,T)  =  z{— - 1)  +  1  or  \ 

=  J(^)3-6(^)2+6(^)-’  pokts  along 

°  0  0  V  AR  n« 

MEMBRANE 

RESULTANT 

/  TRESCA 

/  SQUARE 

n.  =  -1  \ 

DISPLACEMENT 

.  1  (p—  T)  , 

w  (x,  t)  =  —  — —  t2  x 

2  uo 

VELOCITY 

(P-1) 

W  =  -  T  X 

uo 

ACCELERATION 

(p— 1) 

W  =  — -  X 

uo 

TIME  _ 

To 

N/A 

DISPLACEMENT 

AT  REST 

HAS  NOT  COME  TO  REST  YET 
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TABLE  E-2.  SUMMARY,  LONG  SHELLS,  HIGH  LOADING.  2 


TYPE 

SHELL  TYPE 

PRESSURE  LOADING  TYPE 

LONG  C2  >  6 

HIGH  LOAD  p>1+  — 

C2 

CONDITIONS 

u  2-  6 

°  tfjp-l) 

OstsI 

U0<  X  «  1 

MOMENT 

RESULTANT 

mx  =  1 

> 

POINT  B 
,  ON 

MEMBRANE 

RESULTANT 

II 

&■ 

c 

>  TRESCA 
'  SQUARE 

) 

DISPLACEMENT 

W  SB  —  lp-1)  T2 

2 

VELOCITY 

W  =  (p-t)  T 

ACCELERATION 

w  =  p-1 

TIMEt 

ro 

N/A 

DISPLACEMENT 

AT  REST 

HAS  NOT  COME  TO  REST  YET 
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TABLE  E-3.  SUMMARY,  LONG  SHELLS,  HIGH  LOADING,  3 


TYPE 

SHELL  TYPE 

AB  PRESSURE  LOADING  TYPE 

LONG  Cz>  6 

HIGH  p  >  1  +— rr  (p#  —  > 

C  2 

CONDITIONS 

o$x<y  I^T^T, 

/  6  /( p-i)3  L  c2  2  ) 

U°  VC2(p-1)  Tl  P  V  <2p  —  3)  i  6  1  I 

0,  defined  in  TABLE  E4,  ATTACHMENT  1 

MOMENT 

RESULTANT 

mx  =  -1 

MEMBRANE 

RESULTANT 

»V  =  o 

DISPLACEMENT 

1  (P-U  _ 

ui  —  —  A 

2 

VELOCITY 

w  =  0 

ACCELERATION 

w  =  0 

TIMEt 

To 

N/A 

DISPLACEMENT 

AT  REST 

HAS  NOT  COME  TO  REST  YET 
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TABLE  E-4.  SUMMARY,  LONG  SHELLS,  HIGH  LOADING,  4 


TYPE 

SHELL  TYPE 

AB  PRESSURE  LOADING  TYPE 

LONG  C2>  6 

SEE 

HIGH  LOAD  p>l+JL  (p#— ) 

C2  2 

CONDITIONS 

ATTACHMENT  1. 

y<x<u  1  <t  <t, 

0>-  «  2-I2p~3,<P-t|’.  C  ) 

C2  (_  (p-1)3  J 

'-VteSi-T-.i 

m,n(l,  -f).  p>2 

MOMENT 

RESULTANT 

m  |x,t)-  l3  f  4x3+6  (u+yix2  12  uyx 

*  (u-  y|  L 

-  (y+u)(y2-4yu+u2)J 

2  V  3  „  3 

— c  <0,  <  1.  if  p  < 

MEMBRANE 

RESULTANT 

■V--1 

DISPLACEMENT 

SEE  ATTACHMENT  2. 

VELOCITY 

W  -  Ip-T)  Q 

ACCELERATION 

w  =  (V)  1  3  fu  +  y  2x1  1 

\C  /  (u-y3)L  J 

TIME  , 

To 

N/A 

DISPLACEMENT 

AT  REST 

HAS  NOT  COME  TO  REST  YET 
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TABLE  E-4.  ATTACHMENT  1 


) 

) 
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TABLE  E-5.  SUMMARY,  LONG  SHELLS,  HIGH  LOADING,  5 


TYPE 

SHELL  TYPE 

AB  PRESSURE  LOADING  TYPE 

LONG  C2>  6 

HIGH  p>i+-|t  (p*|) 

CONDITIONS 

u<  x<  1 

MOMENT 

RESULTANT 

mx=l 

MEMBRANE 

RESULTANT 

*— 

l 

ii 

9 

e 

DISPLACEMENT 

SEE  ATTACHMENT  1. 

VELOCITY 

W=p-T 

ACCELERATION 

w=-1 

TIME,. 

T0 

N/A 

DISPLACEMENT 
AT  REST 

HAS  NOT  COME  TO  REST  YET 
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TABLE  E-5.  ATTACHMENT  1 

u<x^1  u=u*  atT=x*,  0  =  0* 

w(x,t)  =  —  (p-1)  +  P  [t-t*]  [t2-t,2]  + 

2  2 
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TABLE  E-6.  SUMMARY,  LONG  SHELLS,  HIGH  LOADING,  6 


SHELL  TYPE 


LONG  C2  >  6 


PRESSURE  LOADING  TYPE 


.  6  3 

HIGH  p >  1  +  —  Pft  — 


,  2vy  2 vT 

±V3  ,  {6i  +  C  1  <uo“  C  * 

V,  =  uD  - 9 .  +  —r  lofle  - — - — 

2C  2v^3T  2vT, 

(0 - )  (u  +  ) 

7  C  °  C 


Cy/p~^\ 


r  2^i  /  \ 

Tx/T  «>th  (A)  +  — —  C(1  -  u  J 

^  r  k»  VV 

C  1+  - coth(A) 

L  Cuo  J  (i 


MOMENT 

RESULTANT 


MEMBRANE 

RESULTANT 


DISPLACEMENT 


VELOCITY 


ACCELERATION 


t,  =p- 


(p-DJ  r2_  » 

i-3)  L  6 


0  <  x<  y  y 


,  v  1  (P-1) 

w(x,  r)  = - x 

2 


’  ’  H-f] 


DISPLACEMENT 

AT 

REST 
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TABLE  E-7.  SUMMARY,  LONG  SHELLS,  HIGH  LOADING,  7 


PRESSURE  LOADING  TYPE 


SHELL  TYPE 


LONG  C2  >  6 


2v/3"  2  vT 

V3  .  (0i  +  C  *  (uo“  C  * 

y,  =  -e  + - log*  - 

'  2C  2  3-  2  T 

' lu°  *  — 1 


HIGH  p>1  +  — 


r,  <  r  <  r 

I  O 


r  2V31  ( 

coth  (A)  +  -  C\1  1 

2V3  r _ cu0j  a  = 


VT 

y  ^  x<1  y  <  1  — _ 

C 


r  2V3  -i '  v  J 

««h(A)J  r  c*!!* 

L  -o  J  ; -  2 - d2\ 

1  |  (2p  —  3)  L  6  J| 

,  6  1  (1_vl)  /  \ 

y  « /TT\H)  +  y’ 


02  M2J 


y(  t)  =  1 


-A  1  1  1  ~yi 

c2  (’ -v,)2  ¥7. 


MOMENT  AND 

MEMBRANE 

RESULTANT 


=  1  + - - _  [-4x3  +  6  (1  +  y)  x2  -  12  yx  -  2  (1  -  3y>] 

(1  -V>3  „  _  , 


\  =  “I 


DISPLACEMENT 


w(y,  r)  = 

1  (p  - 1) 

2  U 

o 

1 

(ci-Bi  t) 

T  l09e 

Bi 

(C1-bi"i) 

(1  -x)  1  1 

B1  (C1-B1T)”  (°l  ~BlTl) 


(SEE  NEXT  PAGE 
FOR  MORE 
INFORMATION.) 


ACCELERATION 


To=T1  + 


2(x  -  y)  1 

!1  -  y  )3  ( 1  -  y)2 


[— • $ 


DISPLACEMENT 
AT  REST 


w(yQ  ,  t0) 


when  y  =  yQ  =  1 - the  velocity  becomes  zero. 

C 


(SEE  CONTINUATION 
PAGE.) 
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